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Abstract 

We introduce a new class of Green-Naghdi type models for the propagation of internal waves 
between two (1 + l)-dimensional layers of homogeneous, immiscible, ideal, incompressible, 
irrotational fluids, vertically delimited by a flat bottom and a rigid lid. These models are 
tailored to improve the frequency dispersion of the original bi-layer Green-Naghdi model, and 
in particular to manage high-frequency Kelvin-Helmholtz instabilities, while maintaining its 
precision in the sense of consistency. Our models preserve the Hamiltonian structure, symmetry 
groups and conserved quantities of the original model. We provide a rigorous justification 
of a class of our models thanks to consistency, well-posedness and stability results. These 
results apply in particular to the original Green-Naghdi model as well as to the Saint-Venant 
(hydrostatic shallow-water) system with surface tension. 


1 Introduction 

1.1 Motivation 

This work is dedicated to the study of a bi-fluidic system which consists in two layers of homoge¬ 
neous, immiscible, ideal and incompressible fluids under only the external influence of gravity. Such 
a configuration is commonly used in oceanography, where variations of temperature and salinity 
induce a density stratification; see m and references therein. 

A striking property of the above setting, in contrast with the water-wave case (namely only 
one layer of homogeneous fluid with a free surface) is that the Cauchy problem for the governing 
equations is ill-posed outside of the analytic framework when surface tension is neglected [ 23 ., 2Bi iHQj ■ 
This ill-posedness is caused by the formation of high-frequency [i.e. small wavelength) Kelvin- 
Helmholtz instabilities which are triggered by any non-trivial velocity shear. Recently, Lannes m 
showed that a small amount of surface tension is sufficient to durably regularize the high-frequency 
component of the flow, while the main low-frequency component remains mostly unaffected. 

This result explains why, occasionally, surface tension may be harmlessly neglected in asymptotic 
models, that is simplified models constructed from smallness assumptions on physical properties 
of the flow. This is typically expected to be the case for models derived in the so-called shallow- 
water regime, which implies that the main component of the flow is located at low frequencies; 
and in particular for the two-layer extension of the classical Green-Naghdi model introduced by 
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Miyata [35] [39] , Mal’tseva |36| and Choi and Camassa m- However, as noticed in [29] a linear 
stability analysis indicates that the bi-fluidic Green-Naghdi system actually overestimates Kelvin- 
Helmholtz instabilities, in the sense that the threshold on the velocity shear above which high- 
frequency instabilities are triggered is always smaller for the model than for the full Euler system. 

Many attempts have been made in order to “regularize” the Green-Naghdi model, that is propos¬ 
ing new models with formally the same precision as the original model, but which are not subject 
to high-frequency Kelvin-Helmholtz instabilities, even without surface tension na M 13 ill EE- 
The strategies adopted in these works rely on change of unknowns and/or Benjamin-Bona-Mahony 
type tricks; see EH Section 5.2] for a thorough presentation of such methods in the water-wave 
setting. In this work, we present a new class of modified Green-Naghdi systems obtained through 
a different, somewhat simpler mean. We find numerous advantages in our method: 

• The original Green-Naghdi model is only lightly modified, and the physical interpretations of 
variables and identities of the original model are preserved. 

• The method is quite flexible. It allows in particular to construct models which completely 
suppress large-frequency Kelvin-Helmholtz instabilities; or a model which conforms perfectly 
with the linear stability analysis of the full Euler system. 

• The rich structure of the original Green-Naghdi system (Hamiltonian formulation, groups of 
symmetry, conserved quantities) is maintained. This is generally not the case when change of 
unknowns or BBM tricks are involved; see discussion in EH [53]. 

Our models may be viewed as Green-Naghdi systems with improved frequency dispersion. In 
particular, one of our models shares with full dispersion models such as in m in 03 mi the 
property that their dispersion relation is the same as the one of the full Euler system. Notice 
however that, consistently with the derivation of the original Green-Naghdi model and contrarily 
to the aforementioned ones, our models do not rely on an assumption of small amplitude: we only 
assume that the wavelength is large compared with the depth of the two layers (see the consistency 
result given in Proposition |5.7[ ). To the best of our knowledge, there does not exist any other model 
in the literature with such properties, even when only one layer of fluid is involved. A noteworthy 
feature of our models is that, by construction, they involve non-local operators (Fourier multipliers). 
Such operators are common in deep-water models or small-amplitude models, such as the Benjamin- 
Ono or Whitham equations for instance, but appear to be original in the shallow-water setting. 

Our new class of models is obtained by regularizing the original Green-Naghdi one and not from 
a direct derivation from the full Euler system. As it is made clear in Section |2.1| below, we could 
also derive our models from Hamilton’s principle, but once again our approximate Hamiltonian 
functional consists only in a harmless perturbation of the original Green-Naghdi one. It would 
be interesting to derive our models directly from an expansion of the so-called Dirichlet-Neumann 
operators, and investigate whether the full dispersion model is more precise than the original one 
when only small nonlinearities are involved. 

We also acknowledge that the use of surface tension in view of modeling miscible fluids of 
different densities, such as water and brine, is arguable. We view surface tension as a simple and 
effective way of regularizing the flow. Let us point out that the construction of our new models 
does not rely on any particular structure of the surface tension component, thus could be applied 
with any additive regularizing term. Our systems may be more evidently applicable to two-layer 
systems of genuinely immiscible, homogeneous and inviscid fluids. We believe they may also be of 
interest in the situation of one layer of homogeneous fluid, with or without surface tension. 

For the sake of simplicity, our study is restricted to the setting of a flat bottom, rigid lid and one¬ 
dimensional horizontal variable. The construction of our models, however, may straightforwardly 
be extended to the two-dimensional case. We also expect that our strategy can be favorably applied 
to more general configurations (non-trivial topography, free surface, multi-layer, etc.) 
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In the present work, we motivate our models through the study of Kelvin-Helmholtz instabilities, 
by linearizing the systems around solutions with constant shear. This formal study is supported 
by numerical simulations, which demonstrate how the predictions of the modified Green-Naghdi 
models may vary dramatically depending on their large-frequency dispersion properties, and the 
significant influence of small surface tension. We also provide a rigorous analysis for a class of our 
models by proving the well-posedness of the Cauchy problem as well as consistency and stability 
results, which together offer the full justification of our asymptotic models, in the sense described 
in m ■ This includes the original Green-Naghdi model as well as the Saint-Venant (hydrostatic 
shallow-water) system with surface tension, for one or two layers of fluid; all these results being 
original as far as we know. 


1.2 The full Euler system 

For the sake of completeness and in order to fix the notations, we briefly recall the governing 
equations of a two-layer flow in our configuration, that we call full Euler system. We let the 
interested reader refer to mmm for more details. 

The setting consists in two layers (infinite in the horizontal variable, vertically delimited by a 
flat rigid lid and a flat bottom) of immiscible, homogeneous, ideal, incompressible and irrotational 
fluid under only the external influence of gravity. We assume that the interface between the two 
layers is given as the graph of a function, £(t,x), so that the domain of the two fluids at time t is 

^1 = f {0,y), C (t,x) <z< di}, f4 = f -d 2 < z< C (t,x)}. 

Here and thereafter, the subscript i = 1 (resp. i = 2) always refer to the upper (resp. lower) layer. 
The fluids being irrotational, we consider the velocity potentials in each layer, that we denote fa. 
Finally, Pi denotes the pressure inside each layer. 

Let a be the maximum amplitude of the deformation of the interface. We denote by A a 
characteristic horizontal length, say the wavelength of the interface. The celerity of infinitely long 
and small internal waves is given by 


co = \ 9- 


(P2 ^ Pi)did 2 


P2di + pid 2 


where d\ (resp. d 2 ) is the depth of the upper (resp. lower) layer and p\ (resp. p 2 ) its mass density. 
g denotes the acceleration of gravity. Consequently, we introduce the dimensionless variables 


~ def Z 

z = "T' 

d\ 


~ def X 

X = A’ 


7 def Co 

t - Jt, 


the dimensionless unknowns 

C (t,x) = f ^ t,x \ fa(t,x,z) = f ~^—fa(t,x,z), Pi(t, x, z) c = —— ^2 .Pj(i, x, z) [i = 1,2), 
a aXco clp2Cq 

as well as the following dimensionless parameters 


def Pi 

7 = 5 

P2 


def G. def d-t _ def d\ 

t = Te "“w 5 = dr 


Bo drf g(p2 - Pi)^ 2 


The last parameter is the Bond number, and measures the ratio of gravity forces over capillary 
forces (er is the surface tension coefficient). After applying the above scaling, but withdrawing the 
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tildes for the sake of readability, the system may be written as 


( 1 . 1 ) 


h-d^i + d%(f>i = 0 
d z (j> i = 0 
d z (t> 2 = 0 

( d t C = + ne 2 \d x C,\ 2 d n <j)i = + ne 2 \d x C\ 2 d n (/)2 

a,<h + = 


in SI* (* = 1, 2), 
on {(x,z),z = 1}, 
on {(x,z),z = -6 -1 }, 
on {(x,z),z = eC(t,x)}, 
in (i = 1, 2), 
on {(x,z),z = eC(t,x)}, 


where k(d x Q d = -d x ( ^/ 1+ |^ c|2 d xC) , lP(t,x)J d = hm (p(t, x, <(i, x) + x) - P(t, x, <(i, x) - ;*)) , 

v £ iZ = f (■ y/jxdx,d z ) T and (yi + | 2= « c = -^{d x (){d x <t>i) | a=e? + (d z <j>i) | 2=eC . 

We may conveniently rewrite the above system as two evolution equations, thanks to the use of 

def 

Dirichlet-Neumann operators. In order to do so, we define ip(t,x) = <f>i(t,x,e<^(t,x)), and 


G^ip = G^[eQip d = ^/l + [i\ed x C\ 2 (d n <f>i) \ z=e( = -/J-e(d x Q(d x (j) 1 ) | 2=eC + (d^i) | i=eC 
H»’ s 1> = H»’ s [eC]il> = (02) L=. c = <h(t> 


where <f> i and (f> 2 are uniquely defined (up to an additive constant for (f> 2 ) as the solutions of the 
Laplace problems implied by (1.11. The full system then becomes 


( 1 . 2 ) 


dtC - -G^ = 0, 


d t (d x H^ s il) - ~/d x ip^ + (7 + S)d x ( + ^d x (jdatH^ip ] 2 - 7 |< 9 X ^>| 


= ^d x u^ s - 


7 + 5 d x (k{ty/jid x Q) 


Bo 


c ~VT l 


where 


(iG'0 + eO,C )(0,H>"‘i,)) 2 - 7 (iG^ + e O,0O^)) 2 
M |eC ’*' =-wTTPuP)- 


1.3 Our new class of modified Green-Naghdi models 

Let us now present our new class of models, which aim at describing the flow in the situation where 
the shallowness parameter is small: 

H <C 1. 

Our systems use layer-averaged horizontal velocities as unknowns, that is defining 

If 1 1 r e t 

Ui(t,x) = —r-—r / d x (/) 1 (t,x,z)dz, u 2 (t,x) = —r / d x (/> 2 {t,x,z)dz. 

hi{t,x) J ec h 2 {t,x) J_ s -i 

Here and thereafter, hi = 1 — (resp. h 2 = <5 _1 + e£) always denotes the depth of the upper 
(resp. lower) layer and is assumed to be positive. One benefit of such a choice of unknowns is the 
exact identities (in contrast with 0{^i 2 ) approximations) due to mass conservation (see e.g. |221 
Proposition 3 and (23)]): 

(1.3) d t hi + ed x {hiUi) = 0 (i = l,2). 

These identities are then supplemented with the following 0(fj, 2 ) approximations: 

(1.4) d t (ui + fJ.Qi[hi]u. j) + + ^<9z(l^| 2 ) = ned x (TZi[hi,Ui\) - —' ^—d x Pi, 
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where P 2 - P 1 = -^ dx (^==L— dx c) and 


Q^[hi]ui = Ui}, ^[huui] d = \{hid x ^Ui ) 2 + Ih- 1 ^ d x F?{hld x F>*Ui}, 

where F(* F i(y/J!D) (i = 1,2) is a Fourier multiplier: 


f)v = f 

The choice of the Fourier multipliers is free, although natural properties for our purpose include 
Fj(0) = 1 and F'(0) = 0 (so that Ff —Id is formally of size 0(fj,)), F, (fc) = F, (|fc|) and 0 < F, : < 1. A 
class of Fourier multipliers for which our rigorous results hold is precised in Definition |5.1[ thereafter, 
and we present three relevant examples below. 

• F'^(^fjiD) = 1 yields the classical two-layer Green-Naghdi model introduced in [3Bi; SUITS], 
supplemented with the surface tension component. 

. F^(^JID) = (1 + ^IHI 2 )- 1 / 2 , with 9 t > 0 is an operator of order —1, and sj~fld x F9 is a 
bounded operator in L 2 , uniformly with respect to /i > 0. As a consequence, this choice yields 
a well-posed system for sufficiently small and regular data, even in absence of surface tension. 


-, with convention <5i = 1,(52 = <5. 


(y/JiD) Y 5 l l ^/Jl\D\tanh(S i 1 ^fJl\D\) S i V|£>| 2 ’ 

The modified Green-Naghdi system with this choice conforms perfectly with the full Euler 
system, as far as the linear stability analysis (see Section [3| is concerned. In particular, its 
dispersion relation is the same as the one of the full Euler system. One may thus hope for an 
improved precision when only weak nonlinearities (e <C 1 in addition to /r <C 1) are involved. 
More precisely, we expect that our model is precise, in the sense of consistency, at the order 
0(/r 2 e) instead of 0(/x 2 ) for the original Green-Naghdi system. 


In ( |1.4| ), Pi (or equivalently P 2 ) plays the role of a Lagrange multiplier enforcing the constraint 
resulting from (1.3) and hi + h 2 = 1 + (5 _1 = Cstt, namely 


(1.5) 


d x (hiui) + d x (h 2 u 2 ) = 0. 


Notice now that in the one-dimensional setting, supplementing the above equation with the condi¬ 
tions at infinity Iq — > 0 (|ru| — > 00), one obtains the identity 


( 1 . 6 ) 


hiui + h 2 u 2 = 0. 


In this situation, one may equivalently rewrite (1.3)-(1.4) as two scalar evolution equations: 
dtC + d x w = 0, 


(1.7) 


d t 


hi + 7/12 
h \ h 2 


w + /jQ[e(\v?J + (j + S)d x C + 


hl ~ lh2 'w\ 2 


h 2 y 

= ned x (H[e(,w}) + ^< 9 2 ( 


V^i+^ 2 |a,CI : 


;d x C 


def 
W = 


^ 1 ^’ 2 (u 2 - 7U1) = -hiUi = h 2 u 2 


where, by ( 1.61 and since hi,h 2 > 0 
( 1 . 8 ) 
and 

(1.9) Q F [eC]w d = Q 2 [h 2 ]{h 2 1 w)-')Q\[hi]{-hi 1 w), Tl F [eC„w\ A = Tl 2 [h 2 ,h 2 1 w\-'yTl F i[hi,-hi 1 w\. 


hi + "ih 2 
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Our models are all precise in the sense of consistency at the same order as the original Green- 
Naghdi model, namely 0{^l 2 ). Furthermore, we fully justify system (1.71 for a class of Fourier 


multipliers including all of the three examples above, in the sense that for sufficiently small and 
regular initial data, the model admits a unique strong solution which is proved to be close (when 
/i <C 1) to the solution of the full Euler system with corresponding initial data. However, only for 
the latter two examples is the smallness assumption ensuring the stability of the flow, by preventing 
the appearance of Kelvin-Helmholtz instabilities, consistent with the one of the full Euler system. 
The precise statements of these results are displayed in Section [5] 


1.4 Outline of the paper 

Some elementary properties of our models are studied in Section [2] More precisely, we show 
that all of our models enjoy a Hamiltonian structure, symmetry groups and conserved quantities, 
consistently with the already known properties of the original Green-Naghdi model and the full 
Euler system. 

In Section [3] we recall the linear analysis of Kelvin-Helmholtz instabilities for the full Euler 
system, and extend the study to our models. In particular, we recover that the classical Green- 
Naghdi model overestimates Kelvin-Helmholtz instabilities, whereas our modified model with the 
choice Fj = F‘ mp recovers perfectly the behavior of the full Euler system. 

Section [4] is dedicated to numerical illustrations of this phenomenon. We give two examples 
(with and without surface tension) where the original, improved and regularized Green-Naghdi 
models predict very different behavior. Roughly speaking, the flows are very similar as long as no 
instabilities are present, but the threshold above which Kelvin-Helmholtz instabilities are triggered 
varies dramatically from one model to the other. 

Section [5] is dedicated to the rigorous justification of our models, for a class of admissible 
Fourier multipliers F, . We first prove that the Cauchy problem for system (1.71 with sufficiently 


regular initial data is well-posed under some hyperbolicity conditions. Roughly speaking, we show 
that provided a dimensionless parameter which depends on the large frequency behavior of F, is 
sufficiently small, our system is well-posed for sufficiently regular and bounded initial data, on a 
time interval uniform with respect to compact sets of parameters; see Theorem |5.3| for details. 
We then supplement this result with consistency (Proposition |5.7| ) and stability (Proposition 5.8 1 
results, which together offer the full justification of our models (Proposition |5.9| ). 

Finally, we present in Section [6] some improved results in the limiting case /i = 0, that is on the 
so-called Saint-Venant, or shallow-water system (with surface tension). 

Appendix [A] is dedicated to the detailed presentation of our functional setting and some nota¬ 
tions; and Appendix [B] provides the proof of our main result, Theorem |5.3| 


2 Hamiltonian structure, group of symmetries and conserved 
quantities 

It is known since the seminal work of Zakharov m that the full Euler system (with one layer) 
admits a Hamiltonian structure. This Hamiltonian structure has been extended to the two-layer case 
in 0(35] 28]. We show in Section 2.1 that our models - both under the form fHijl-fTdj) and ( fTfl ) 
also admit a Hamiltonian structure, so that our models could be derived from Hamilton’s principle 
on an approximate Lagrangian. Such an approach has been worked out in [37] [I2J I13| (see also [Mj) 
in the one-layer case, in m in the regime of small-amplitude long waves or small steepness, in in 
the two-layer case with free surface, and lacks for existing regularized Green-Naghdi systems in the 


literature mM®m- We then enumerate the group of symmetries of the system (Section |2.2| ) 
that originates from the full Euler system (see (61), and deduce the related conserved quantities 
(Section [23| . 
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2.1 Hamiltonian formulation 

Let us define the functional 


( 2 - 1 ) 

'H F (h 1 ,h 2 ,u 1 ,u 2 ) = \ f -4y~^ ~—{h\-^h\-C)+-(h 1 +h 2 -H)P 1 + 2<yl +J 1 ^ (y/l + fi\d x h 2 \ 2 -l) 
2 J R e 2 1 - 7 e fie 2 Bo 

+ 7^il w i| 2 + h 2 \u 2 \ 2 + fi^hi(hid x F^ui) 2 + fi^h 2 {h 2 d x F^u 2 ) 2 

with H = 1 + <5 _1 (the total height) and C = cq + c\h\ + c 2 h 2 (Casimir invariants) where Co, ci, c 2 
are chosen such that the integral is well-defined. It is convenient for our purpose to define 


def , _ 

Wi = hiUi 


and 


Vi =* A F [hi)wi =* -j-Wi - ^h i 1 d x F^{hfd x F^{h i 1 w i }}. 


Under the non-cavitation assumption hi > 0, A F [hi] is a symmetric, coercive thus positive definite 
operator, and therefore we may equivalently write the Hamiltonian functional as 


/ H F {h 1 ,h 2 ,Vi,v 2 ) = \ ( -^j^-{h 2 ~"fh 2 -C)+-(h 1 +h 2 -H)P 1 + 2 ^ 1 ^ (i/l + fi\d x h 2 \ 2 -l) 
2 J R e 2 1 - 7 e fie 2 Bo 

+ 7 ui* 4 i[/ii] -1 ui + v 2 A F {h 2 \~ 1 v 2 . 


We deduce the functional derivatives 

sn F , Fr , lH1 , _ sn F 

- — ='yA 1 [h 1 \ Vi="fh 1 ui, —— 

OV 1 0V 2 

^2 + c i) + ~ Fl + 2 l^ 1 ! 2 “ 


— — h 2 U 2 , 


Shi 

sn F 

Sh 2 


1 7f{,, . 1 7 + 5 „ 

(h 2 + c 2 ) -|—Pi -„ „ <9, 


d x h 2 


e 2 1 — 7 


e 2 Bo V^i + A ,|a x/l2 |2 


+ ~ /r7^2[^2, U 2 ]. 


One can now observe that the system (1.31-( 1.4) enjoys the non-canonical symplectic form 


d t U + eJ^— = 0 with U = f 
oU 


Vi 

h 2 


\ 

( 7 1 d x 

\ 

and J = f 

T l d x 

dx 

) 

\ 

d x J 


Let us now indicate how the Hamiltonian formulation of our system can be identified with the 
well-known similar formulation of the full Euler system for internal waves, which we recall below. 
We base our discussion on the two-equations system ( |1.7| ) as it allows a direct comparison with 
earlier works on the bi-fluidic full Euler system, such as |19| . We introduce the functional 

( 2 . 2 ) m,v) = l [ (7 + S)C 2 + (Vl + ^ 2 M 2 - 1) - «Si(7& + Qi)- l Q 2 v, 

/ j R fie rso 

where Qi = Gi[e£] are such that d x (Gid x tp) = Gip, and Gi are Dirichlet-to-Neumann operators: 

• G\ : H> L( d n q i) |_ =ef , where (j> is the unique solution to 

fid 2 (j) + d 2 (j) = 0 in </>(x,e£) = <p and d z <f>(x, 1) = 0. 


• G 2 : H> i(9 n ^) |_ =ef , where <j) is the unique solution to 

fid 2 (f) + d 2 (j) = 0 in Oj, </>(£,£<() — and d z (j){x, — S~ F ) = 0. 
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The operators Gi are well-defined if h 2 {eC) > h 0 > 0 (see [321 Chapter 3]), and consequently 

Gid x <p as well, thanks to the identity |321 Prop. 3.35] 

Gi[e(]ip = -d x (hiUi) (* = 1,2). 

Now, recall the construction of the full Euler in Section |1.2| and in particular the definitions of 
(f>i(t,x,z),(f> 2 (t,x,z) and ip(t,x) = x, e£(t, x)). Denoting 


def 

v(t,x) = d x ((j) 2 (t,x,eC(t,x)) ~'y^i(t,x,eC(t,x))), 

we deduce the identities v = G 2 l G\d x Y — id-xY = ^ 2 _1 (l?i + iG^dx'tp. 

With these definitions, one can check that the full Euler system (1.21 can be written as 

( 8H' 


d t C = ~d x 


V 5v 


d t v = -d x 


fsn 


We now can view the Hamiltonian functional of our system as a 0{fi 2 ) approximation of the 
Hamiltonian functional of the full Euler system. Indeed, one has (see e.g. (32 Prop. 3.37] when 
F 1 - = 1, but the general case adds only a 0(fi 2 ) perturbation): 

QiHT 1 = [hi]+0(n 2 ), Ai[hi] : w ^ ^-w ~ {hfd^ {h^w}}. 

Notice that, using the definitions of Section [131 

A F [e(\w d = jA[[h!]w + A F 2 [h 2 }w = j-w + - "/fi^-h^ 1 d x F^{hld x F^{h]; 1 w}} - fi^-h^ 1 d x Fl^ {hld x F^{h 2 1 w}} 

tl 2 O O 


+ jfiQ F 1 [hi](h 1 1 w) + fiQ F [h 2 )(h 2 1 w) = hl+ ' yh2 w + fiQ F [eQ 

n\il2 


W. 


In particular 7l F [e((] is a symmetric, coercive positive definite operator, if hi(eC)-, h 2 {eC) > ho > 0. 
Plugging the above (truncated) approximation in the full Euler’s Hamiltonian functional \2.2\ yields 

(2.3) n F (C,v) d ^l [ (7 + 5)( 2 + (v/1 + ^ 2 |9,C| 2 - 1) + vAAQ-A 

z j R fie no 

which corresponds to the previously defined Hamiltonian functional ( |2.1[ ) when hi + h 2 = H = 1 + (5 _1 
Indeed, we may introduce 


(2.4) 


w= f 7l F [eC] l v and m = f (— l) l w/hi (i = 1,2), 


so that 

A{C,v)= \ f (7 + *)C 2 + 2i V^W± + ^ 2 |5,C| 2 - 1) + wA F [eC]w 
z j R fie no 

= \ [ (7 + ^)C 2 + W l + M£ 2 I^CI 2 - 1) + "lhi\u\\ 2 + h 2 \u2\ 2 

z j R fie Jjo 

+ A^hi{hid x F^uiY + fi — h 2 {h 2 d x ^2 U 2 ) 2 . 

Computing the functional derivatives, we get 


sn F 

Sv 

sn F 

sc 


7l F [eC] 1 v = w and 


(7 + S)C 


l±Aa ( (3*0 

Bo x V(l + /xe 2 |^C | 2 ) 1/2 


eh\ — ^/h 2 
2 h\h 2 2 


wl 


and one recognizes the Hamiltonian structure of system 0>: 

'sn F ' 


dtC = ~dx 


5v 


d t v = -d x 



fieTZ F [e C,w], 
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2.2 Symmetry groups 

Based on the work of [ 6 ], one may list symmetry groups of our systems. Most of the symmetry 
groups of the full Euler system have no equivalent for the Green-Naghdi model, because they involve 
variations on the vertical variable, which is not accessible anymore. The most physical symmetries, 
however, remain. We list them below. 


Assume (C,u,,Pi) (with i = 1,2) satisfies (1.3)-(1.4). Then for any k £ M, (£■* ,uf , P*) also 
satisfies (1.3)-(1.4), where 


i. Horizontal translation 

(C,u?,Pn(t,x) = f (c(t,x- x),Ui{t,x- x),Pi{t,x - *:)) 

ii. Time translation 

(C,P?){t,x) d = (C(t- K,x),Ui{t- 3t,x),Pi(t- x,x)) 

iii. Variation of base-level for potential pressure 

d = (c(t,x),Ui(t,x),Pi{t,x) + x) 

iv. Horizontal Galilean boost 

(C, ufi p n (*> x ) = f (C {t, X - xt) , Ui (t, x - Kt) + xr, Pi (t, x - xt)^j 


It is interesting to notice that when working with formulation (1.7). i.,ii.,iii. induce symmetry 


groups as well (although iii. is trivial), but not iv.. Indeed, because the Galilean boost breaks the 
conditions Uj —> 0 at infinity, the identity w = —h\Ui = I 12 U 2 cannot be deduced anymore from the 
constraint coming from the rigid-lid assumption. Thus the set of admissible solutions to ( |1.7| ) 
is too restrictive to allow Galilean invariance. 

2.3 Conserved quantities 

The first obviously conserved quantity of our systems, given by (1.31, is the (excess of) mass: 

- 77 ^ = 0 , Z(t) d = f C(t,x)dx. 

at JR 


Equations (1.41 yield other conserved quantities: the “horizontal velocity mass” 

-V» = 0, Vi = f [ Ui + nQ^[eC]uidx (i = 1,2). 
c Jr 


d_. 
d t' 


Choi and Camassa m observed a similar conservation law of the original Green-Naghdi model, 
and related this result to the irrotationality assumption of the full Euler system. Indeed, one has 


H « / d x (<j)i(t,x,eC(t,x)))dx, 
Jr 


where the velocity potentials <pi (i = 1,2) have been defined in Section |l.2| and the approximation 
is meant with precision d(/i 2 ). Thus one has by construction = C(/r“), and it turns out that 
this approximately conserved quantity is actually exactly conserved by the Green-Naghdi flow (see 
also [23). Of course, the linear combination (recall (|1.8[) and (1.9 l) 


V = f V 2 - 7V1 = f \ + ^ w + pQ F [eC] w da: 
Jr hi h-2 
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is a conserved quantity of system (1.71. The above conserved quantities may be interpreted as 


Casimir invariants of the Hamiltonian system. 

After long but straightforward manipulations, one may check 


l M = 


h\d x P\ + h,2d x P2 — [h\Pi 4- ^ 2 ^ 2 ]!^ 


+ OO 
OO ’ 


M d = 


jh±ui + h,2U2dx. 


The total horizontal momentum is in general not conserved. This somewhat unintuitive result is a 
consequence of the rigid-lid assumption and has been thoroughly studied in mm- 
One has the conservation of total energy: 


A^ f -0, - j^n + 5)C + (y/l + ne 2 \d x C\ 2 - 1 ) 

r y 1 

+ "/hi\ui \ 2 + h 2 \u 2 \ 2 + fi-^hi(hid x FiUi ) 2 + ti-h 2 (h2d x F^u 2 ) dx. 

The conservation of energy may be deduced from the Hamiltonian structure of the system and more 
specifically from the invariance of the Hamiltonian with respect to time translation; see e.g. m- 
Similarly, the invariance with respect to horizontal space translation yields the conservation of the 
horizontal impulse: 

~T,' I = 0, Z(C,v)= f (v dx. 

Jr 

This conserved quantity of the bi-fluidic Green-Naghdi model seems to have been unnoticed until 
now. In the one-layer case free surface, that is when 7 = 0, the horizontal impulse is related to the 
momentum (which is conserved in this situation) through the horizontal velocity mass: 

< 5 _1 V + el = / h.2V dx = / h 2 U2 + Q2[ e C]^2 da; = / h 2 U2 dx = M. 

J R J R J K 


In this situation, the symmetry with respect to Galilean boost yields an additional conserved 
quantity, which is the counterpart of the “horizontal coordinate of mass centroid times mass” for 
the full Euler system as defined in [B], namely 


C(t) d ^ f 


C(t,x)dx, with C(t, x) = (x — th, 2 V or, equivalently, C(t,x) = (x — tw. 


The conservation of C can be deduced from the above, as (1.31 yields 

—C = [ xdtC — w — td t w d.T = [ —xd x w — w da: — t— [ wdx = —t—M. = 0 . 

dt Jr Jr dt Jr dt 

Remark 2.1 (Traveling waves). The Hamiltonian structure and conserved quantities of our system 
provide a natural ground for the study of traveling wave solutions. For instance, one easily checks 


that critical points of the functional H f — yield traveling wave solutions to system (1.7). 

However, as soon as the surface tension component or non-local operators are present, explicit 
formula such as the one provided in JM \EB for the original Green-Naghdi model seem hopeless. 
We thus postpone to a further study the analysis of traveling wave solutions to our system. 


3 Kelvin-Helmholtz instabilities 

In this section, we formally investigate the conditions for the appearance of Kelvin-Helmholtz 
instabilities for the full Euler system as well as for our Green-Naghdi models, based on the linear 
analysis of infinitely small disturbances of steady states. We say that the system suffers from Kelvin- 
Helmholtz instabilities when arbitrarily high frequency unstable modes are predicted by this linear 
theory. 
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3.1 The full Euler system 


We linearize (1.1 ) around the constant shear solution: (£ = 0 + x(, (f>i = UiX+xcfi i, </> 2 = u 2 x + x(f> 2 ) 
where ^ and u 2 are constants, and x«l. Notice that by (1.51, one has necessarily u 1 + S~ 1 u 2 = 0; 

and therefore Uj = and u 2 = ypj, where v *= f u 2 — 7 Mi- When withdrawing 0(x 2 ) terms, one 
obtains the following linear system (see e.g. [33]): 


(3.1) 


dt( + c(D)d x ( + b(D)d x v = 0, 
d t v + a(D)d x ( + c(D)d x v = 0 


where v = d x ((4> 2 - 7^i) |, =eC ), and 

^ ^tanh(^//I|fc|) - 7 tanh(fyfyW 1 |A:|) ev 
c \k) — 


tanh( v /^|fc|) + 7 tanh( % //.i< 5 _ 1 |fc|) 7 + 5’ 


b(k) d =l f 1 


tanh(y^|fc|) tanh( v //l(5 1 \k\) 
y/Jl\k\ tanh(y/p|fc|) + 7 tanh(y / /r 5 _ 1 |/c|) 


and 


%/p|fc|7 


(5 + 1) 2 


a ( k ) (7 + <^)( 1 + go ) tanh(^l|A:|)+ 7 tanh( v ^ 5 - 1 |fc|) (<5 + 7 ) 2|e - r 


Since b(k) > 0, the mode with wavenumber k is stable (namely the planewave solutions e^ fcx u ±( k '> t '> 
satisfy u±(k) £ K) if and only if a{k) > 0. For small values of k, this yields the necessary condition 

2| , 2 S(S + l ) 2 , . 

76 tel 7jT^ <7 + i 

For large values of k, one approximates tanh(fy7t|fc|) + 7 tanh(y / /Z 5 _ 1 |fc|) « 1 + 7 , and deduce 

/ n\i ~ / , r\ 7 2 MBo (5 + 1) 4 4| |4 

mm {«(*)}« (7 + S) - 4(1+7)2(7 + ii)s E l»l ■ 

The full Euler system is therefore stable for each wavenumber provided 

T|v | 2 ^ 7(1 + \Jn Bo)e 2 |u | 2 is sufficiently small. 

3.2 Our class of Green-Naghdi systems 


When linearizing (1.71 around the constant shear solution, U (£, w) T f = (0 + x£,w + xw) T , 
where w is constant, one obtains the following system: 


(3.2) 


dtC + 9 x u> = 0, 


b ( D)d t w - c F (D)d t ( + a F (D)d x C + c F (D)d x w = 0, 


with 




35 


and 


-F/nJ'f/ | c\ 21 |2/c3 , \ 21 i 2 ^1 ^2 1 2 + l \^1 I~ i 1 12 . 7 + ^i,i 2 

a (k) = (7 + 5) - e \w\ (5 + 7 ) - /re \w\ --- \k\ + 


where (with a slight abuse of notations) Ff* = F i{JJik). 
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The stability criterion is more easily seen when rewriting system (3.2 1 with unknown 0 


Indeed, one obtains in that case 
(3.3) 
with 


35 

dtC + c F (D)d x ( + b F (D)d x v = 0, 
d t v + a F (D)d x ( + c F (D)d x v = 0 


c F {k) = ew ——pr 


S a -y | .. |Fg| 3 -7|F5T |,|2 

+ ^ 3(7+5) \ K \ 




1 + f-l 


35(7+5) 


b F (k) = 


1 

7+<5 


|FgP+7<|Ff| a MU2 


1 + n 


35(7+5) 


and 


/m = ( 7 + 5) + Ifcl 2 - Icuf 7((5 + 1)2 ^ + ^l fc l 2 l F 2l 2 )( 1 + 5M|fc| 2 |F^| 2 ) 

{ ) h ’ bo 11 M 5( 7 +5) i+ ^' 2 s yi[ F ; 12 \ k ? 


35(7+5) 

Analogously with the full Euler system, the mode with wavenumber k is st able if and only if 
a F (fc) > 0. Let us quickly discuss the three examples introduced in Section 


1.3 


• In the case of the original Green-Naghdi system, F\ d (y/~jID) = 1, the condition to ensure that 
all modes are stable is 

YgnM 2 = f 7(1 + MBo)e 2 |w | 2 is sufficiently small. 

This is more stringent than the corresponding condition on the full Euler system in the 
oceanographic context, where one expects )(Bo> 1, as noticed in [251155] , 


If F; mp (^D) = 


(with convention 5i = 1,52 =5), 


y 5- 1 VMl^|tanh(5- 1 VM|!5|) ^Vl^l 2 

then the linearized system ( |3.3| ) is exactly ( |3.1[ ) (recall that by (1.5), w = ; A ; ju): 


a F (k) = a(k) ; b F (k) = b(k) ; c F (fc) = c(k). 


In particular, the stability criterion of this Green-Naghdi model corresponds to the one of 
the full Euler system. As previously mentioned, this also shows that the model has the same 
dispersion relation as the full Euler system, as this corresponds to setting w = 0. 

• In the case F‘ eg (,/Z ID) = — , - , one remarks that 

1 v /1 + MIDP’ 


(7 + 5) > |ew| 


! 7(5 + 1 ) 2 (5 2 + _L )(i + J—) 

5(7 + 5) 1 30 2 n 30/’ 


is a sufficient condition to ensure that all modes are stable, and does not require the presence 
of surf ace tension. A natural choice is 6 i = with convention S i = 1 ,82 — 8 , for the Taylor 
expansion of the dispersion relation around /i = 0 to fit the one of the full Euler system, at 
augmented order 0(/x 3 ), instead of the 0(/x 2 ) precision of the original Green-Naghdi system. 

lr This change of unknown is not without signification. It consists in writing the system with the “original” variables 
of the full Euler system: v = d x {(<l >2 — 7</>i) | z=e £ ) ? or more precisely G(fi 2 ) approximations of these variables. It 
is interesting to note that the nonlinear hyperbolicity condition in Theorem |5.3| matches the naive sufficient condition 
for stability a F (k) > 0 coming from ( |3.2| ) but not the sharp condition a F (fc) > 0 (in particular when 7 —> 0); see also 
Remark |5.5| 
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In Figure [l] we plot the instability curves corresponding to a F (fc) for the three above examples. 
More precisely, for fixed k £ R, we plot the value of e 2 |re | 2 above which a F (k) > 0, and thus 
instabilities are triggered. One clearly sees a great discrepancy for large wavenumbers. In particular 
the minimum of the curve, which corresponds to the domain where all wavenumbers are stable, not 
only varies for each model, but also is obtained at different values of k. 



Figure 1 : Instability curves of the modified Green-Naghdi models with F, = 1 (original), F, = F‘ eg 
(regularized) and F^ = F‘ mp (improved). The last one coincides with the full Euler system counter¬ 
part. The dimensionless parameters are 7 = 0.95, <5 = 0.5, e = 0.5, /i = 0.1, Bo ^ 1 = 5 x 10 -4 . 


4 Numerical illustrations 

We numerically compute several of our Green-Naghdi systems, with and without surface tension, 
in order to observe how the different frequency dispersion may affect the appearance of Kelvin- 
Helmholtz instabilities. 

As in Figure [l] we focus on the three aforementioned examples: F^ = 1 (original), F^ = F‘ eg 
with 6 i = (regularized) and F, = F‘ mp (improved). Values for the dimensionless parameters 

are as in Figure [l] above. The initial data is C(0,a:) = — e -4 ^” and w(0,x) = 0. 

Figures [2] and [3] represent the predicted flow at time t = 2 and t = 3, in the situation with 
surface tension. Figure [4] represents the predicted flow at time t = 2 in the situation without 
surface tension. Each time, the left panel plots the flux, w(t, x) (or rather 1 + w for the sake of 
readability) as well as the interface deformation, ((t,x ); while the right panel plots the spatial 
Fourier transform of the interface deformation, £(t, k). The dashed line represents the initial data, 
and the three colored lines the predictions of each model. 

Discussion In the situation with surface tension, we see that at time t = 2 (Figure [2]), the predic¬ 
tions of the three models are similar. Only the original model shows small but clear discrepancy, and 
most importantly early signs of instabilities. This situation is clearer when looking at the Fourier 
transform, right panel. We see the existence of a strong large frequency component which has grown 
from machine precision noise. As expected, modes with higher wavenumbers grow faster. Notice 
the regularized model also exhibits a non-trivial (although very small) high-frequency component. 
This component is however stable with respect to time: it is not produced by Kelvin-Helmholtz 
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Figure 2: Prediction of the Green-Naghdi models, with surface tension, at time t = 2. 




Figure 3: Prediction of the Green-Naghdi models, with surface tension, at time t = 3. 




Figure 4: Prediction of the Green-Naghdi models, without surface tension, at time t = 2. 
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instabilities, but rather by numerical errors. It does not appear when surface tension is absent 
(Figure [I]) . 

At later time t = 3 (Figure [3]), the Kelvin-Helmholtz instabilities have completely destroyed the 
flow of the original model. For such irregular data, spectral methods are completely inappropri¬ 
ate and we do not claim that our numerical scheme computes any relevant approximate solution. 
Meanwhile, the flows predicted by the regularized and improved models remain smooth and are very 
similar. When running the numerical simulation for much larger time, our computations indicate 
that the flow of the regularized and improved models remains smooth for all time. 

When surface tension is neglected from the models, we see (Figure [4]) that at time t = 2, Kelvin- 
Helmholtz instabilities have already destroyed the flow for the original model, and the improved 
model shows early signs of instabilities in its high-frequency component. The flow predicted by the 
regularized model, however, remains smooth and is very similar to the flow with surface tension. 


Numerical scheme Let us now briefly present our numerical scheme. It is very natural in our 
context to use spectral methods [49] as for the space discretization, since Fourier multipliers are 
dealt similarly as regular differential operators. Such methods yield an exponential accuracy with 
respect to the spatial mesh size for smooth data. In our simulations, we used 2 9 = 512 equally 
distributed points (with periodic boundary conditions) on x £ [—4; 4]. As for the time evolution, 
we use the Matlab solver ode45, which is based on the fourth and fifth order Runge-Kutta-Merson 
method [49, with a relative tolerance of 10~ 10 and absolute tolerance of 10~ 12 . It is convenient 
to solve the system written in terms of £ and v = w + fiQ F [e(]w; see (2.41, although this 

requires to solve at each time step w as a function of £ and v. 

In Table [l] we display the numerical variations, between time t = 1 and initial time t = 0, of the 
conserved quantities (discussed in Section 2.3) as a very rough mean to appreciate the precision of 
the numerical scheme. One sees that the agreement is excellent, except when the horizontal impulse 
is concerned. The latter shows a great sensibility to the presence of large frequency components, 
indicating that such component, as expected, affects the precision of the numerical scheme. It is 
remarkable that the other conserved quantities do not suffer from such a loss of precision. 



With surface tension 

Without surface tension 


original 

regularized 

improved 

original 

regularized 

improved 

Mass Z 

1.5321 icr 14 

1.5654 10“ 14 

1.5876 10" 14 

1.299 10" 14 

1.5321 10 -14 

1.4988 10 -14 

Velocity V 

-6.3135 10- 16 

-7.3647 10- 17 

9.1148 10“ 16 

-1.4962e-15 lO^ 15 

3.7576 lO" 16 

5.3084 10 -16 

Impulse X 

5.4362 10“ 10 

6.6559 10" 17 

-2.1199 10- 17 

-1.3741 10“ 2 

-4.4256 lO” 17 

2.2034 10~ 14 

Energy £ 

-3.5282 10~ 12 

-7.3064 10~ 13 

-4.8171 10- 12 

-3.3628 10“ 12 

-5.1278 lO’ 12 

-4.7436 lO’ 12 


Table 1: Difference between conserved quantities at time t = 2 and time t = 0. 


5 Full justification 


This section is dedicated to the proof of the main results of this work, namely the rigorous justifi¬ 
cation of the class of modified Green-Naghdi systems introduced in Section [TT3| i.e. 


dtC + d x w = 0, 

(5T) dt (h^ w + ^ [<]w 


+ (7 + fi)d x C + f ^{jhl 


-iK 




= + ^ d K vw\o*c\* dxC )' 
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where 


Q F [e{]w = -\h^d x ^ 2 {hld x ^{h^w}} - ^dxF^hld^ih^w}}, 
n F [e(,w] = ^wh^dxF^hldxF^h^w}} - ^wh^dxf^hldxF^h^w}} 

+ \(h 2 d x F%{h 2 1 w}) 2 - ^{h.dxF^w}) 2 . 

Here and thereafter, we always denote h\ = h\{eC,) = 1 — e( and /12 = ^(eCf) = £> _1 + e(- Let us 
also recall that F/ (* = 1,2) denotes a Fourier multiplier: 

Ff = F. UVi^D) i.e. Fj7(0 = Fi(V^)/(0- 

In order to allow for the functional analysis detailed in Section |B.1| we restrict ourselves to admis¬ 
sible Fourier multipliers F/, in the following sense. 

Definition 5.1. The operator Ff = F (i = 1,2) is admissible if it satisfies: 

i. F,; : R — > R + is even and positive; 

ii. Fi is of twice differentiable, F,;(0) = 1, F'(0) = 0 and sup fcgR |F"(fc)| < C f < 00 ; 

in. k 1 —y |fc|F,(fc) is sub-additive, namely for any k,l £ R, \k + Z|F,;(fc + l) < |fc|F,(fc) + |Z|F,(Z). 

In that case, one can define appropriate pairs Kp. £ R + and c, £ [0,1] such that 


(5.2) 


Vfc £ R, Fi(fc) < K Fi \k\~*. 


Proposition 5.2. The three examples provided in Section 1.3 namely F‘ d , Fj eg , F) mp are admissible, 
and satisfy (5.2 ) with (respectively) <; = 0,1,1/2. 

Proof. A sufficient condition for F/ to be admissible is to satisfy, in addition to i. and ii., 
iii’. k 1 —^ fcFj(fc) is non-decreasing on R + and k 1 —> F,-(fc) is non-increasing on R + . 


Indeed, for k + l > k > l > 0, that F, (fc) is non-increasing on R + yields 0 < F^k + l) < Fj (k) < F, (Z), 
and therefore (k + l)F,(k + l) < fcF^fc) + IFifk) < kF^k) + IF^l). This shows k £ R + 1 —> kFi(k) = 
\k\Fi(k) is sub-additive. Since F, is even, one shows in the same way that k £ R _ 1 —>■ \k\Fi(k) is 
sub-additive. Finally, for k < 0 < l, that F, is even and k £ R + fcF, is non-decreasing yields 
|k + l\Fi(k + l) = \k + Z|F,(|fc + 1 1) < (|fc| + |Z|)Fj(|fe| + |Z|), and the sub-additivity in R + allows to 
conclude. 

That property iii’ holds is immediate for the first two examples, only the last one requires 
clarifications. One easily checks that k £ R + 1 —>■ fcF) mp (fc) is non-decreasing, so that we focus on 
the proof that k 1 —>■ F/ np (/c) is non-increasing for k £ R + . To this aim, it suffices to show that 


Vfc > 0, f(k) =* — k 2 — ^fcsinh(2/c) + cosh(2fc) — 1 < 0. 


This follows from /(0) = /'(0) = /"(0) = 0 and 


f"'(k ) = 2sinh(2fc) — 4fccosh(2fc) = 2 cosh(2/c)(tanh(2/c) — 2k) <0 (fc > 0). 


Thus property iii’ holds, and k |fc|F‘ rap (fc) is sub-additive. 


□: 
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For reasons explained in Appendix [B] our energy space involves both space and time derivatives 
of the unknowns. For U = (£,«;), we define 


N 


E"(U) d = |C 


A' 0 


+ \ w \ 


E n (U) = E°(d a U) = |C 
M=o 


\w\ 


where N always denotes an integer and a a multi-index. The functional setting and in particular 
the definitions of functional spaces A'^-i and are given in Appendix [a] 

In addition to 7 ,/r, e, <5, Bo -1 > 0, it is convenient to introduce the dimensionless parameters 


(5.3) 


T f d = e 2 (l + ( 7 ^ Fl + A'f 2 )(pBo) 1 -^) < 00 


where <;, Kf 1 , Kf 2 are specified in Definition 5.1 (5.2); and 


def 


m = max{e, 7 , 5, 6 ,/r,Bo } < 00 . 


Theorem 5.3 (Well-posedness). Let U° '= (£ °,w°) T £ x Y^f with N > 4, satisfying 

h] d = 1 - eC° > h 0 > 0 , h° 2 d = <5 _1 + eC° > h 0 > 0 ; 

( 7 + 5) - e 2 max {((h°)~ 3 + 7 (/i?)- 3 )kT} > k 0 > 0 . 


One can define K = C^m, h 0 , k 0 1 , e|£| ff3 ) such that if Tf|u>° 


z i <K 

1 


-1 


there exists T > 0 and 


a unique U ^ (£,w) T £ C^([0, T); X^ 0 -i x solution to (5.1) and U | t=0 = U°. Moreover, 
there exists C 0 = C(m, h^ 1 , k^ 1 , K, E N (U 0 )) such that 


T - 1 < C 0 x 


( e + T^>\ FV +T F |u; 0 |^) 


and sup E N (U) < C 0 X E N (U°). 

te[o,T) 


The proof of this result is postponed to Appendix [B] Let us for now discuss a few implications. 
Remark 5.4 (Initial data). Since our functional spaces involve time derivatives, it is not a priori 


clear how to define C u 


lx* 


and ke 0 


1 Cl 


As it is manifest from the proof the definition of 


(d a U°) | (=0 for sufficiently regular C l °(x),w°(x) is given by system (5.1) itself. More precisely, for 
a = (0,02); then the definition is clear. We then define ( d a U° ) | t=0 for a = (01,0:2) with. a\ > 0 
by finite induction on oi, through the identities obtained from differentiating 1 oi| — 1 times 

with respect to time. These identities given in Lemma \BA\ and are uniquely solved by Lemma \B. b\ 

Remark 5.5 (Domain of hyperbolicity and time of existence). Hypotheses on the initial data 
ensure that the flow lies in the “domain of hyperbolicity’’ of the system; see Lemma B.5\ They 
may be seen as the nonlinear version of the stability criterion presented in Section \3.S\ as they 
provide sufficient conditions for Kelvin-Helmholtz instabilities not to appear. However, remark that 
our “Kelvin-Helmholtz instability parameter”, T F , is not multiplied by 7, in contrast with T and 
Tqn In Section [5| as well as the nonlinear criterion on the full Euler system given by Lannes m 
(5.1)]. The latter results imply that the large-frequency Kelvin-Helmholtz instabilities disappear in 
the limit 7 —► 0, so that surface tension is not necessary for the well-posedness of the system when 
7 = 0 . We do not recover such property with our rigorous analysis, although numerical simulations 
indicate that our models are well-posed when 7 = 0 and Bo = 00, as long as the non-vanishing 
depth condition is satisfied. 

A second setback is that the time of existence involves Ty 2 |u >°|_ 1 , and not only TF|ru °|„ 2 • 
In practice, this means that when I F < 1, and in particular when T F < e C 1, then the time of 
existence of our result is significantly smaller than the one in \32l Theorem 6 ]. 
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However, let us note that our conclusions, in particular with the choice F; = F‘ mp where <; = 1/2, 
are in complete agreement with aforementioned results in the oceanographic setting of internal waves, 
where one expects large values of e and 7 « 1. 

We believe that the above limitations originate from the choice of unknowns used when quasilin- 
earizing the equation. This was quickly discussed in footnote^in Section\3.2\as for the occurrence of 


7. The restriction on the time of existence originates from estimates (B.8) and (B.9) in Lemma B .4 


and more precisely the lack of an analogue of m Lemma 7] thanks to which “good unknowns” were 
constructed. We show in Section^how the techniques used in this work, applied to the Saint-Venant 
system (that is setting p = 0) written with different unknowns, yields sharp results. 


Remark 5.6 (Regularized systems). In the case g = 1, one sees that Theorem \5.3\ does not depend 
on Bo (through Tf). In particular, the results hold true even when surface tension is neglected, 
i.e. Bo -1 = 0, and we recover in that case the “quasilinear timescale” T -1 < e. Our strategy relying 
on the use of space-time energy is not needed in that case, as classical energy methods can be applied 
to prove the well-posedness for initial data in Sobolev spaces: (<( 0 ,u> 0 ) T £ H 8 x H s , s > 3/2. 

We show now that the above well-posedness analysis can be supplemented with consistency and 
stability results, which together provide the full justification of our models, (5.11. Such a program 
was completed for similar models in the one-layer setting in OHIO ( see I32j for a detailed account). 
In the two-layer case with rigid upper lid, the consistency of many models were derived in [5]. The 
consistency result below builds upon m and m- 


Proposition 5.7 (Consistency). Let U (</ ,ip) T be a solution of the full Euler system (1.21 such 

that such that there exists C$,T > 0 with 


Icl 


L°°([0,T);HZ +5 ) 


l^cll 


\d x %b\\ w + \\d t d x 'ib\ 


y ^ Co, 
L°°([0,T);ff x +2 ) 


for given s > to + 1/2, to > 1/2- Moreover, assume that there exists ho > 0 such that 

(5.4) \/(t,x) € [0,T) x R, hi(t,x) = 1 — ef(t,x) > ho > 0, h 2 (t,x ) = 5 -1 +e((t,x) > h 0 > 0. 


Define w by d x w = — ^G^\eC,\4> = —d t Q. Then (/’,m) T satisfies exactly the first equation of (5.11, 
and the second up to a remainder, r, bounded as 

IMIl oo ([o,t);.h' < ') — h x C(m, h 0 , Cq, C f). 

Proof. The Proposition has been stated and proved, in the case of the original Green-Naghdi system, 
F‘ d = 1, in (20) Proposition 2.4]. By triangular inequality, there only remains to estimate 


fQ = f h\dt(Q pd [eC\w - Q F [eCM 
Notice that by Definition |5.1| one has 


Hi 


and r-jz = pe\d x [eC,w] - H F [eC,w}) 


Hi 


|Fj( v / /ffc) — l| < -C F /r|fc| 2 


It follows, since H 8 is an algebra for s > 1/2 and by Lemma B.3 that 


r e<ESl^K rla 4 F ?- Id )K^ F .WM})l ff . + ?I^K 7l ^K^( F r-ld)R- 1 t«}}) 


HV 


i=l 

,2 


< p C F C(m, h 0 i\dtC\ L aa^ 0T y H s+4y\d tW \ Loo ^ 0 T).H^+ i )y 


r iz S p L f C(m,h 0 ) |C| i oo([ 0i T);ffi +5 )> I w L°°([o,t) ; ^ +6 ))' 


Similarly, one has 
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Now, recall that by definition, d x w = —dt( £ L°°([0, T); Hf +4 ), and w £ L°°([0, T); L 2 ) by 
identity and the uniform control of up, see e.g. [22] Proposition 4], Thus we control 

w £ I/°°([0, T); H x +5 ), and similarly d t w £ L°°([0, T); H x +i ), and the Proposition follows. P 

Proposition 5.8 (Stability). Let N > 4 and U, = ((i,Wi) T £ L°°([0, T); X^ o _ x x Ypf) solution 
to (5.11 with remainder terms (0,ri) T . Assume Q satisfy ( |5.4[ ) ,(B.27),( BT28] ) with ho,ko,K > 0. 
Set 2 < n < N — 1 and assume that d a ri £ ^([(LT); (Ip),)*) for any |a| < n. 

Then there exists 0 < T* < T such that for all t £ [0, T*) 

E n (U 1 - U 2 ) 1/2 < C 0 E n (U 1 | t=0 - U 2 U„ ) 1 / 2 e xt + C 0 f f n (t')dt'. 

Jo 


with 


A = Co x + Tf|| t/ii| 


L“([0,T];Z f A) 


), fn(t) = l 5 “ ri “ S “ r2 l(Y^)*’ 


|of| <n 


and C 0 = C(m,h 0 ,k 0 \K,\\Ui \\ Loo ^ 0 T] , x4 iXY * t y II^ 2 IIl°°([o,t];X 4 iX y F ^))- Moreover, one has 
[T ) < C 0 (e + T F ||^i|| i oo([ 0iT ]; Z i (i ) + ^f|| u,1 IIl°°([o,t];z;L))' 


Proof. By Lemma B.4 for any |a| < n < N — 1, = f (d a (i,d a Wi) T satisfies (B.30) with 

remainder terms fj a ^ r^ + d a r i £ L 1 ([0,T); (Y^)*), and 

\r[ a) ~r ( 2 a) \. Y0 < \d a n - d a r 2 \ ( y 0 V +C 0 x (e + T^ /2 | Wl | +X F | Wl || 1 ) x E™ (lh - U*) 1 ' 2 , 


with C' 0 = C(m, ho,E N {U 1 ), E N {U 2 )). By Lemma 


B.10 


one has 


E°(u[ a) - C/ 2 (a °) 1/2 < Co E°{u[ a) \ t=0 -U^ a) \ t= 0 ) 1 / 2 e xt + C 0 f / ( “ ) (^ , )e A(^ ■ t ' ) d^ , , 

■A) 


with Cq, A as in the statement and 


f ia \t) = | f\ a) - 


(“) A«) | 




+ e l C72 l(w^°°) 2 l^ 1 ~ 1,2 1 


Since n > 2 and ./V > 4, one can restrict T* as in the statement and augment Co if necessary so 
that the estimate holds. □ 

The following is now a straightforward consequence of Theorem |5.3| Propositions |5.7| and |5.8| 


Proposition 5.9 (Full justification). Let U° = ((°,w°) T £ Xft a -i x Y F ^ with N sufficiently large, 
and satisfying the hypotheses of Theorem 5.3 Define ip° with d x w° = and assume 

that (<C°,'0 O ) T satisfies the hypotheses of Theorem 5 in JETf- Then there exists C,T > 0 such that 


There exists a unique solution U = to the full Euler system (1.21, defined on [0,7*] 

and with initial data (<^°, ■0°) T (provided by Theorem 5 in JSTj); 


There exists a unique solution Up = (^F)'*i , f) T to our modified Gree n-N aghdi model (5.1), 
defined on [0,T] and with initial data (<(°,w°) T (provided by Theorem 5.3); 

— ^G^[eC,\ip = —dtC, one has for any t £ [0,T], 


|(C> W ) (Cf, w F)|| ZjOO ([ 0it ]; A -0 _ iX yo 1 ) — C v t - 


• Defining d x w = 
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6 The Saint-Venant system 


The Saint-Venant system (with surface tension) is obtained from our Green-Naghdi models (5.1 1 
by setting fj, = 0. The results of Section [5] thus apply as a particular case. However, it is possible 
to obtain sharper results by considering the system obtained after the following change of variable 

— def h\-\-~*h 2 — — 

v = = u 2 - 7«i: 


( 6 . 1 ) 


dtC + d x (H{e Qv) = 0 , 

d t v + (7 + 6 )d x ( + f3 x (fT(eC)M 2 ) - = 0 , 


where we denote H(e <) = h '^l 2 ■ 

In the following, we quickly review the steps of the method developed in Appendix [B] providing 
results without proof. 

The analogue of Lemma |B.4| is the following: 


Lemma 6.1. Let U = (C,^) G -^b - 1 x with N > 2, solution to (6.11 and satisfying 


( 6 . 2 ) 


hi(e() = 1 - e( > h 0 > 0, h 2 (e() =5 + e( > h 0 > 0. 


For any a = (a.\,a 2 ) such that |a| < N, denote U ^ = f (d a (,d a v) T andv^ ^ (d a ei v,d a e 2 v) T 
(if a.j = 0, then d a ~ ej v = 0 by convention). Then U^ satisfies: 

dtC ia) + d x (H(eOv {a) ) + eH'(e()vd x (^ + 3 x 0 a [eC]h<“> = r [ a) , 

d t vW + cUsvK,™]C (a) + eH’{eC,)vd x v^ = r { 2 a) , 

withf> a [e(}v {dl l d = i, 2 } a j- ff, ( e 0 ( e5e3 'C)(^“ _ei «) and 

cisv [<,«!]• d = ((7 + S)+ e —H"{eCfi\v\ 2 ^ 

and eh] = (r[ a ^ [e£, eh], [e£, eh]) T G -^b 0 -i x L 2 satisfies 


.1+1 d V 
Bo 


r-(“) 


[ e C> eh] 


K„-i xL2 


< e C(m, h 0 1 , |C 


\x e 


, \v 




r) X (|C| 


X' 


I H* 


One has immediately the following analogue of Lemma |B.5| 

Lemma 6.2. Let (£, h) T G L°° x L°° he such that satisfies (6.2) with h 0 > 0, and 


(6.3) 


(hi + h 2 ) 2 




hi 2 > ko > 0. 


Then there exists Kq,K\ = C(m,h 0 1 , A : 0 1 , eld .„) such that 


V/, 5 GL 2 , 


V/eA-“ 0 _„ 


V/Gl 2 


| (&sv[eC? ex j ]/, g) . p | — |h|y° 5 

^ Bo -12 Bo -1 Bo -1 

\(H(eQf,g) L2 \ < A' 1 |/| i2 |gf| i2 , 

1 2 

<a sv [ee,eh]/,/) (xl _ i)t >—l/l^, 

(if(eC)/,/) i2 >^|/| 2 i2 . 
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A priori energy estimates are obtained by adding the L 2 inner product of the first equation 
with asvC^ Q \ and the one of second one with H{eQv^ a ' 1 + and following the proof of 

Applying the strategy of Section |B.5| one then obtains the following 


Lemmata |B.8| and B.10 
analogue of Theorem n 


Theorem 6.3. Let N > 2 andU 0 *= f (£°,u°) T £ X^ o _ 1 xH N , satisfying ( |6.2[ ) , ( 6.31 with ho, ko > 0. 
Then there exists T > 0 and a unique solution U ^ (£,u) T 
ing (6.11 . Moreover, there exists Co = C(m, h$ , 


M^°l; 


: C«([0 ) r);X ^_ 1 X H n ) satisfy- 
xH n) such that 


T” 1 < C, 


< Cq x e and 


sup 

te[o,r) 


c 




I H> 


< Co 




lx* 


I® 0 1 


H N 


Remark 6.4. Theorem 6.3 is valid uniformly with respect to the parameter Bo” 1 , and the result 
holds in particular in the case without surface tension: Bo” 1 = 0. This case is however straightfor¬ 
ward as the Saint- Venant system is then a quasilinear system, and the result was stated in particular 
in ]2&f . Assumption (6.3) corresponds exactly to the hyperbolicity condition provided therein. 

0: Kelvin-Helmholtz instabilities 


Notice that (6.31 is automatically satisfied in the limit 7 


disappear in the water-wave setting, and ( 6 . 2 ) suffices to ensure the stability of the flow. 


Remark 6.5. While completing our work, we have been informed of related results covering the 


water-wave setting (y = 0,<5 = \) of our Theorem 6.3 Saut, Wang and Xu in w generalize the 
result to Boussinesq systems in dimension d = 2, while Chiron and Benzoni-Gavage m treat more 
general Euler-Korteweg systems. 


A Notations and functional setting 

The notation a < b means that a < Co b, where Co is a nonnegative constant whose exact expres¬ 
sion is of no importance. We denote by C(Ai, \ 2 ,...) a nonnegative constant depending on the 
parameters Ai, A2,... and whose dependence on the A j is always assumed to be nondecreasing. 

In this paper, we sometimes work with norms involving derivatives in both space and time 
variables. We find it convenient to use the following sometimes non-standard notations. 

• For 1 < p < 00, we denote L p = L p = L P (R) the standard Lebesgue spaces associated with 
the norm 

I fW = ^ \f(x)\ p dx^j < 00. 

The real inner product of any functions fi and /2 in the Hilbert space L 2 (R) is denoted by 

(fij2) L 2= [ fi(x)f 2 {x)dx. 

J K 

The space L°° = Lff = L°°(M) consists of all essentially bounded, Lebesgue-measurable 
functions / with the norm 

I/Loo = esssup 2 , gR |/(a:)| < 00. 

• We acknowledge the fact that only space derivatives are involved by the use of a subscript. 
For k £ N, we denote by W^ , 00 (]R) = {/, s.t. VO < j < k, d 3 x f £ L°°(]R)} endowed with its 
canonical norm. 

For any real constant s £ R, = H*(M.) denotes the Sobolev space of all tempered distri¬ 
butions / with the norm |/|#» = |A S /|^2 < 00, where A is the pseudo-differential operator 
A = (1 - 3 2 ) 1 / 2 . We denote Hff = f) NeN H*. 
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• In absence of subscript, the derivatives are with respect to space and time, and thus apply 
to functions defined on (t, x) £ [0, T ) x R. Thus for N £ N, W N, °° is the space of functions 
endowed with the following norm: 

I/U- = £ |9"/l 

\ot\<N 

where we use the standard multi-index notation: a € N 2 , < 9 ( ai >“ 2 ) = d^d ^ 2 and |a| = ai + a 2 - 
In particular, d ei = c^ 1,0 ) = d t and d e2 = = d x . 

Similarly, H N is the space of functions endowed with 

Ifl 2 = V' ia“f | 2 

\ J \ H N Z^ I J II 2 ■ 

[a|<JV 


We denote H°° = fljveN-f *. 


• Given p, 7 , Bo 1 > 0 and F, (i = 1,2) admissible functions (in the sense of Definition 5.1 1 , 
we define F f ° , W ° o _ 15 Z% as the completion of the Schwartz space, <S(R), for the 

following norms: 


1/1 


.Y° 


def I „ 1 2 1 I „ , 1 2 

= Wv + Bo&f 


def 


L 2 ’ 


/ v° = / r2 +M7 \dx^lf\r2 +M \d x F%f 


I L 2 ’ 


\f\ W o =\f\ L i+il\dxf\ L i, \f\ z o =\f\ L 1 +Vin\W\ L l+M&ff\Li- 

Bo - l -DO FM 

For iVGNwe define consistently with above the norms controlling space and time derivatives: 

2 

yO ? 

- 1 cu 


1 fl 2 d _£ f \ " 1 ao f I 2 

1° J\x° ’ 

|a|<JV 

|/lv« d = E \ d °f\ 

fa|<jV 

1 f 1 de f \ 1 30 f 1 

\J lw* , “ Z^ J \ W° 

\t\ ZN d = E \ da f\ 

l l Z/p M z- J 1 1 


|a|<JV 


|a|<AT 


• Denoting X any of the previously defined functional spaces, we denote by X* its topological 

dual, endowed with the norm |y>| = sup(|y>(/)|, \f\ x < 1); and by (-,-)(x)* the (A* — A) 

duality brackets. 

• For any function u = u(t, x) defined on [0, T) xl with T > 0, and any of the previously 
defined functional spaces, X, we denote L°°([0,T);X) the space of functions such that u(t, •) 
is controlled in A', uniformly for t £ [0 ,T), and denote the associated norm 

IMI roo/m yi = eSS Sup|u(f, •) I < OO. 

II II L ([0 ,T),X) te[0iT) ' 

For k £ N, C k ([0,T); X) denotes the space of X-valued continuous functions on [0,T) with 
continuous derivatives up to the order k. Finally, C^([0, T); A”) is the space of continuous 
functions with values in X, given the weak topology. 


B Proof of the Theorem 15.3 

This section is dedicated to the proof of our main result, Theorem |5.3| The proof relies on energy 
estimates which are also used in the proof of Proposition |5.8| 

Our strategy is similar to the one used for the full Euler system with surface tension by 
Lannes EH ED, and originates from an idea of Rousset and Tzvetkov EH 121 . The main difference 
with respect to the traditional methods for quasilinear systems is that we treat time derivatives in 
the same way as space derivatives. In particular, the main tool of the analysis is the control of a 
space-time energy. The reason for such a strategy is that 







A new class of two-layer Green-Naghdi systems with improved frequency dispersion 


23 


the two unknowns, f and w, are controlled in different functional spaces, one being contin¬ 


uously embedded in the other but to the price of a non-uniform constant (see Lemma B.ll, 
and the inclusion being strict; 

• the most singular term of the system, namely the one which involves the operator of highest 
order, comes from the surface tension component, and couples the two unknowns (it appears 
as an off-diagonal component of the quasilinearized system). 

This is why one cannot use standard energy methods in Sobolev-based functional spaces, as com¬ 
mutator estimates fail to control all coupling terms; see also the discussion in pH] . 

More precisely, our strategy is as follows. In Lemma |B.4| below, we “quasilinearize” the system. 
We differentiate several times the equations with respect to space and time, and extract the leading 
order components. The quasilinear system we consider is the complete system of all the equations 
satisfied by the original unknowns and their space-time derivatives up to sufficiently high order. 
Thus only L 2 -type estimates on the aforementioned linear “block” systems will be required. In 
Section [B.3| we study the operators involved in the block systems, and exhibit sufficient conditions 
for the existence of a coercive symmetrizer of the system in Lemma |B.5| This yields a priori 


energy estimates in Section B.4 Finally, in Section |B.5| we explain how to deduce from these 
energy estimates the well-posedness of the linear block systems (Lemma |B .11 ), and in turn the 
well-posedness of the nonlinear system (Theorem |5.3[). 


B.l Technical results 

In this section, we provide tools (injections and product estimates) similar to the classical ones 
concerning Sobolev spaces, for the functional spaces X^ ) _ 1 , Yp]f, , as defined in Appendix [XJ 


Let us fix p, 7 , Bo 1 > 0 and F, (z = 1,2) admissible functions (in the sense of Definition 5.11. 
In particular there exists K Fi > 0 and £ [0,1] such that 

(B.l) |F i (0| 2 <min{l, J A Fi |£r 2 <}. 

The following standard injections will be frequently used, sometimes without notice: 

(B.2) \f\ Lao < \ f \ L1 < \ f\ H t 0 (t 0 > 1/2); thus \f\ Z N < |/|y~+i and \ f\ Z ” ~ \ f\ H »+*o+i- 

One immediately sees that the space is continuously embedded in V//; the following 

Lemma precises the norm of the inclusion map. 


Lemma B.l. IfFi satisfies (B.ll, then 

I 2 


V/ G X° 


f\ Y o < (1 + (7* Fl + a: F2 )(mB o) w ) |/| 




Proof. The inequality is a simple consequence of Parseval’s identity and Young’s inequality: 

V ^I a )l/l0l a 


li\d x Ffi^fiD)f\l 2 < J \t\ 2 ~ 2 vm 2 d£ < K Fi (v Bo) 

where we used Bo^ |^| 2-2 ^ < sBo+(l — c)^! 2 - 


i-c 


0 


Sobolev spaces Hff and enjoy straightforward product estimates, which are immediately 

N T/T/iV 


extended to and W ^ o _ x : 


\f9 \x° 

Bo - 

and therefore 


< 


ttlw 0 \f\ 


x° 


< 


\X 1 


\x° 


VN > 1, 


\fg\ 


x' 


,-1 

'N 


< I/I 


X' 


\fg\ 


\x r 


w° 


< f 


w° 


\9\w° 


The following Lemma shows that spaces Yp/ and enjoys similar estimates, thanks to the 
sub-additivity property of admissible functions (recall Definition 5.11. 
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Lemma B.2. Let F, 


(i = 1,2) be admissible functions. Then for any 1 < p,q,p,q,r < oo 


satisfying 1 + 1 = 1 + | 

+ 

-Hl'a, 

II 

one 

has 







(B.3) 

\d x mfg} 

\ L r < \f 

1 Lp 

|a 

rFf. 


^ 5 

i/tf' 7 /' 1 

LpI 17 * 1 * J \ L9 

It follows in particular: 










(B.4) 


\fg 

y° 

1 FP 

< 

\g\ 


1/ 

| Y° Z 
*FP 

9 y 1 / V° I 

1 FM 

(B.5) 

\/N> 1, 

\fg 

yN 

< 

|/| 

VlV 

\g\ 

yN ? 
J FP 


(B.6) 


\fg\ 


< 

|/| 


\g 

kv 



Proof. From the sub-additivity, one has ■ v /Ml^l F i(v / A*0 < + Vt |£ - r?|F;(\/ m(£ “ ’?))• 

Thus 


\d x ^{fg}\ r Lr = / mt^)Y\f*g\ r {m= / dC 


<y d , 
< Ut 


d??|^|F - v) 

J Mv\Fi{Vi Ir l)\f\(v)\d\((, - v) + 1C - ??|Fi(\/^(£ - r?))|/|(r?)||g|(^ - 77 ) 

(IWM?l)(0 + (l/M^I)(0 


where we used that \d x F?f\(£) = \i^Fi(y/JI^)f^)\ = |£|F i ( v //Z£)|/(£)| since F i(y/jlQ > 0. Esti¬ 
mate (B.3) follows from Young’s inequality for convolutions. 


Estimate (B.4) is deduced with r = p = q = 2 and p = q = 1, and using 


Estimate (|B.5|) follows from the above result and triangular inequality, 

< 


\fg\ Y »< E c a , 0iN \(d a f)(dP g )\ 


Y° 

1 CLl 


E l«“/lw l^/h 


|a+/3|<JV 


|a|<AT-l,|/3|<iV 


Estimate (B.6) follows from (B.31 with p = p = q = q = r= 1. 


□ 


We now provide Schauder-type estimates in our functional spaces. 

Lemma B.3. Let H £ C°°(— < 5 _1 ,1) and e£ £ L°° such that 

fti(eC) = 1 - < > h 0 > 0, h 2 (e() = d> _1 + e( > h 0 > 0. 

Then, denoting H n ^h 0 = \H\ Cn ^_ s - 1+h 1 _ h ^ and fixing to > 1/2, one has 

• For any s > 0, if C, £ H s and f £ Hf, then one has with n £ N, n > max{s, to}: 


H(eC)f\jja , Hn,h 0 , | e C|^max{s,t 0 }) |/ 




• For any f £ L 1 , one has 


\H{e()f \ L1 < C(h 0 \H 1M , |eC| ff t 0 )|/l L i- 

• For any N £ N, if C £ H t x 0+1+N and f £ Zthen one has 

\H(eQf\ zN <C(ho\H 2 +NM ,\e(\ „t 0 +i+tv)|/| „ N . 

F P x F A 4, 
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• For any N £ N, if ( £ H* 0+1+N and f £ then one has 

|-^(^C)./| ^ C(h v , J^2-\-N,ho 5 | | rt^p+ 1 + N ) | f | y|W - 

Proof. In each case, we decompose H(eQf = H(0)f + (H(e() — H(0)) f = H(0)f + Gh 0 (e()f where 
Gh 0 is such that G^ 0 £ G°°(K), Gh 0 (x) = H(x) — H( 0) for x £ [— J -1 + h 0 , 1 — h 0 ] and Gh 0 {x) = 0 
for x £ R \ [—<5 _1 ,1]. It is clear that, since min{/ii(e£), /z.2(eC)} > ho > 0, one can construct such a 
G ho satisfying additionally: for any n £ N, \G ho \ cn = C(ho 1 , H n>ho ). 

The first estimate is a direct consequence of a classical Schauder-type estimates in Sobolev 
spaces; see e.g. [48]. As for the second, one has 

|^0/Li < \H(W\ Ll + \G^M)*f\ Ll <ff(0)|^ Ll + |G^K)| il |/] il . 

The second estimate now follows from (B.2) and applying the above result: 

|^o( e C)| L i < |^o( e C)| ff *o < C(^0 ,h 0 , | e C| ff *o)- 

Using that Z^,, is an algebra, one has 

\H(eOf\ z o < \H(0)f\ z o + \G h M)f\ z o < m\f\ z o + \G h M)\ z o \f\ z o • 


Since \u\ 


y0 < kl 

Z^r-Ll 


l ff . 0 +i for any u £ LT‘ 0+1 , one deduces the third estimate for N = 0 as above. The 
case N > 1 is obtained by induction, differentiating N times H(eC,)f and applying Leibniz’s rule. 
The last estimate is obtained identically since by Lemma [B. 2 1 


□ 


\H(eOf\ Y o <\H{0)f\ yO + \G ho (eOf\y° < H(0)\f\ yO + \G h M)\ z ° \f\ Y o ■ 

1 -* mi. 1 1 -*■ mi ' 1 i mi. ' ' M mi. ' ' mi. 1 1 x mi. 


The proof is now complete. 


B.2 Quasilinearization of the system 

The following Lemma introduces the quasilinear block systems which are central in our analysis. 


Lemma B.4. Let U = (£, iw) £ Ag o _ 1 x Yf( with N > 4, solution to (5.1 1 and satisfying 


(B.7) 


/ii(eC) = 1 - e£ > h 0 > 0, h 2 (eQ =5 + e( > h 0 > 0. 


For any a = (an, 02) such that \a\ < N, denote U = f {d a (, d a w) T and ^ ( d a ei £, d a e2 () T 
(if ctj = 0, then d a ~ ej ( = 0 by convention). Then one can define £ (Y^)* such that 

dt( {a) + d x w^ = 0, 

b[eC]3 t u;(“) + d x a[e(), ew]C (a) + 9 x a Q [eC]C <a) + c[<, ew\d x wW = r (a) , 


where 


> h i + ihl 


a[eC,ew]« d = ((7 + 6 ) - e 2 " 1 ^ '" 2 |w| 2 ) x • - /ze 2 (di TZl[h 2 ,w\ + 'yd 1 'R.[[h 1 ,w])* 


{hih 2 ) 3 


b[e(]« d = hl + 7/t2 « + fi{Ql[h 2 \+'yQ F 1 [h 1 ])t 
n\ri2 


Bo 


d x 


(1 + fie 2 \d x <(\ 2 ) 3 / 2 


c[eC, eu>]» = f 2 e ^ w x • - iie(fiQ^[h 2 \(w) - 7 dQp[/ii]('u;))« 


{hih 2 y 


- g.e(d 2 U F 2 [h 2 ,w] -'yd 2 U F 1 [h 1 ,w])», 
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with dQf, di72.f, defined in (B. 14l,(B.16),(B.20),(B.211 below; and 


a«K]< w = ^4 I £ 3o 


Bo 


.jfie- 


,(d x d^O(d x O(d x ^ &) ] 

(1 + /-ie 2 \d x (\ 2 ) 5 / 2 


if |a| = N, and 0 otherwise. 


\je{ 1 , 2 } 

Moreover, r= r [e£, etc] satisfies 

(B.8) l r(a) |(y F 0 (l) * < C(m,ho\E N (U)) x E^{Uf/ 2 x (e + T^ /2 |^|^ + T F |w|^), 

and 

(B.9) |r^Ki,e Wl ] -r(“)[eC 2 ,eH | (>•<,)* < C(m,h^ ,E N (U 1 ),E N {U 2 )) x E^ilh - U 2 )^ 2 

x (e + T^ /2 | Wl | , +T F | Wl |" ). 


Proof. The proof simply consists in differentiating a times the Green-Naghdi system (5.11. The 
higher order terms contribute to a, b, c and a, while lower order terms contribute to P a h In the 
following, we explain how the estimates concerning are obtained. 


Contribution from the first order terms, d t ( hl ^J^ 2 w) + fA:( |w| 2 ). 

Applying Leibniz’s rule, one finds 


< B10 » 

with 


h 2 - 7 h 2 , |2 \ , _ 2 


w = -e 


h l+'V h * i ...\2 a , ,^1-7^2 ...a «a„. , _(a) 


ihih 2 ) 3 


\w\ 2 d x d a C + e 


(Mi 2 ) 2 


-wd x d a w + er\ , 


l“l + l / n \ / 2 \ l“l + l 

r[ a) = J2 Y, e"C ,(ft -^>GW(eC) n^Cj t[d^w d = ^ ^ e «C'(/^) r f il/3 i ) ) 


n— 0 f3i,(3j \i—l / \j=l 

where (A, A-) is any n + 2-tuple of multi-index satisfying 


"=° ft,ft 


1 < |AI < ••• < \/3 n \ < M, 0 < |AI < \/3' 2 \ < |a| and — a + (0,1)) 

»=1 1=1 

£ 7 (ft,ft) j s a constant and G^ the n-th derivative of G(X) = hl ^ 1 ]Y^(xf^ = • 

We estimate each of these terms as follows : 

• if |A»| = |a|, then 0 < |AI, • • \0n-i\, Wi\, \P 2 \ < 1, and 


(ft-ft)i 


I L 2 


< iG«">K)i t „ia»-ci L , nw*. 


w\ 


,1=1 


U =1 


otherwise 0 < |/3i|,..., |/3 n |, \/3[\ < |a| — 1, and 




I L 2 


< |G W «)L- ni^'CL. (\d^w\ L ,\o>k\ 


L 2 ) ‘ 
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One has \G^ n \eC)\ LOO < C(m, h^ 1 ) since e(( satisfies ( |B.7[ ); and by Sobolev embedding, 

\d^u\ LOO < | d^u\ Hl < min{|u| x i+i3i, |u| f i+i^i}. 

We deduce immediately that for N > 2 and |a| < N, 

(B.ll) l r i a) |^ ^ C(m,ho 1 ,E N (U)) x E^(U) 1/2 . 


For the second contribution, one deduces from the first equation of (5.11, dt( = —d x u>, that 

m n\ oafl fhi + lh 2 \ h 1 +"/h 2 h\ - ^h\ («) 

(B.12) a o t —,—,- w = —,— - o t o w + ew—- —,— -o x o w + er\ 

\ h\h2 ) h\ri2 h\h2 

where r^' 1 is estimated as above: 


(B.13) 


r-r| i2 < C(m,ho\E N (U)) x £ |a| (t/) 1/2 . 


Contribution from the dispersive terms, pd t {Q ¥ [eQw) — ped x (77 F [e£]w;). 

Define (with a slight abuse of notation with respect to Section |1.3| ) 

(B.14) Q^/i»]u> d = {hfd^ {h^w}}, 

so that Q F [eC]w = Q 2 [h 2 ]w + ^Q\[hi]w. Differentiating a + ei times and using d t ( = —d x w yields 

(B.15) d a d t Qi[hi]w = Q\[hi\d a d t w - (-l)MQf [hi](w)(ed a d x w) + 

where 

(B.16) dQfNM* = \K 2 (d x ^{h*d x Pi{hr'w}}) x . 

- h^d x ^{h 2 d x ^{h- l w} x .} + x .}}; 

and 

4? = E^ ) (^ l O^ F f{(^^ 8 )ft.F{‘{(^v 1 )(^)}} = £c (ft v 3? } » 

Pj Pj 

where is a constant and (f3j) is any 4-tuple of multi-index satisfying 

4 

0 < t/3i|, \fi 2 \, |/33|, |Ai| < |a| and = ol + (1,0), 

3 =1 

We estimate each of these terms by assuming that U = (C,w) T £ S(R) x <S(R), so that for any 
/ £ 5(K), the following identities are immediately valid: 

= {i^hr^d^iid^h^d^ii^hr^id^w)}} , f) L2 
= -[(d^)d x F l f{(d^h- 1 )(d^w)} , d x Fni^hi 1 )/}) L2 - 

The estimates hold as well for U = (£, w) T £ X^ o -i x ^r* (-W — 4) and / £ Yff by density of <S(R) 
in and _X"g x using standard continuity arguments. 
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if |/3i| = |a|, then 0 < |/3 2 |, |/3 3 |, |/3 4 | < 1 and 

Notice that for |/3i| = \a\ > 1 there exists j £ {1, 2} such that ej < /3i and 

I d Pl h.\ 2 = | d^(h~\d^()\ L2 < C(m, V\|eC| wM _ 1>oc ) x e|C| HM . 


L 2 ' 


Now, using several times Lemma B.3 and since d^ 2 h : - = 3hfed^ 2 ( if |/3 2 | = 1 or = /i? 

if I/? 2 1 =0 (and similarly for d^ 3 h~ ), one has 




< CKV 1 , | e C | ) | (^ 3 ^i" 1 ) (d^ 4 w) | z i 

FF 4 

< C(m, /iq x , |eC|^ 4 ) | w |• 

Therefore, since max{4, |a|} < N, 

7 2 ~>(r 3 ^\/) ( r F °j* < C(m, Iiq 1 , |<| HW ) x ey/'f- i n\w\ z2 |C| X M \f\ L 2 - 

Bo 1 

if |/3 2 | = |a|, then 0 < |/3i|, |/3 3 |, |/3 4 | < 1 and 

One has as above d^ 2 hf | J2 < e|^| ff | Q| C'(m, | e C|v^id—i)- By Lemma [fTls] one has 

V'y 2 ~ i v\dxFi{(d Pl K 1 )f }\ L2 < \{d^h£ 1 )f\ Y „^ < |eC| ff 3 )|/| r o i - 

The last term is treated identically and one obtains eventually 

7 2_ V(^\/)(y F °J*| < \|eC| HN ) x e|C| x w_ i Mz 4 m \f\ Y ^. 

• if |/3 4 | = |a|, then 0 < |/?i|, |/3 2 |, |^ 3 | < 1 and 

l 2 - i ^ ) J)(Yo l£ y <7 2 ->|^/ l f|^|^Ff{(^/ l -i)(^ U ;)}| i 2 |a ;c Ff{(^/ l - 1 )/}| 
Reasoning as above and since |/3i| + |/3 2 | + |/3 3 | = 1, one obtains 

\ < CKVMeC| H3 ) x e\<;\ z iJw\ Y w\f\ Y ^. 

• if |/3 3 | = | a |, one obtains as above 

< C(m,hc\ |eC| ffS , |cC| wl a|-i,«) x 


L 2 ‘ 


By Lemma B.l and (5.31, it follows 


7 2 V(r 3 ^\/) ( r F °j* < |eC| ffN ) x Tp /2 |w| zl |C| 


c „, ,.)f |a| , \ J Iyo ’ 

FF 4 Bo -1 FF 4 
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Y° • 
1 FA* 


• otherwise, one has 0 < |/3i|, |/? 2 1 5 1/?31, \Pa\ < H ~ 1. 

If |/3i| < |a| — 2, we integrate by parts and estimate 

7 2 -V<r$ ) ,/>(yo.).| < I|^Fr{(^ 31 )(^ 4 ^)}II^Ff {(a ^ 31 1 )/} I 
If 1/011 = | a | — 1, then we estimate 

7 2 ”V (45 } 5 /) i2 | < M|^ 1 ^- 1 | i2 |^ F f{(a^/ i f)c? a; Ff{(^3 /l - 1 )(^^)}}|^| / | i2 . 

In both cases, we find 

7/>(*&)* < CKV 1 , Kltfiv, My") X e|C| x M_ i l/ly^- 


Plugging these estimates into (B.15I, we proved 
(B.17) 

nd a d t (Q F [e(]w) = n(Q F [h 2 ] + r )Q F 1 [h 1 ])d a d t w - ii(dQ F [h 2 ](w) - 'ydQ F 1 [h 1 ]{w)) (ed a d x w) + r ( 3 a) , 
with 


(B.18) 


r£ a) | (n) * < C(m, ho 1 ,B N (U)) x ^{U ) 1 ' 2 x ( e + T^>| zl ). 


The other contribution is treated similarly. We define 7?. F [e£,u>] = 72-2[^ 2 — 'y7t F [hi,w\ with 

(B.19) n F [KM d ^ f l -wh- 2 d x ^{h*d x ^{h- l w}} + \(h i d x ^{hr 1 w}) 2 , 

(B.20) di TZ F [h i: w]» d = -^wh~ 3 d x F l ^{h 3 d x Ff(h~ 1 w)} x • 

+ wh~ 2 d x F^{h 2 d x F^{h~\u} x • } ^whr 2 d x F?{h%d x F?{hr 2 w}} 

+ (hid x Ff {h^w}) x {[d x F^{h^w}) x • - (hid x F^{h~ 2 w x • })), 

(B.21) d 2 7 Z F [hiM* = l^d^ihld^ih-'w}}) x . + ±whr 2 d x F!?{h%d x F!t{h? x .}} 

+ (M*Ff{/M™}) x x •}). 

It follows 

(B.22) /iei9 Q (9 a; (^ F [eC,w]) = Vt 2 dx((difc 2 [h 2 ,w] + 7 diK F 1 [h 1 ,w])d a (;') 

+ ned x {[d 2 n F 2 [h 2 ,w\ - yd 2 TZ F 1 [h 1 ,w])d a w^ + r[ a) , 

where r M may be estimated similarly as r M above: 

(B.23) < C(m,ho 1 ,E N (U)) x ^(t/) 1 / 2 x (e + T F \w\l^). 
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Contribution from the surface tension term , ^^d 2 ( , 1 cLC). 

1 1 ’ Bo x Vv/i+^ia.ci 2 J 

Let us denote s(d x () = -^=L^=^d x ( and notice ds = ^ 1 +fie 2 \g x ^ 2 ) 3/2 dd x (. It follows 


(B.24) d a d 2 x s = dl 


) -1 ) +4“’. 


(1 + ne 2 \d x C \ 2 ) 3 / 2 U Ux y 1 + ^WI 2 ) 5/2 


with 


N+l / 2 \k 2k-\-l 

r (“) - V _ r W d 2 L f TT d P Jd c 

5 ^ (l + /xe 2 |5 x C| 2 ) fe+3 / 2 rfc ’ xC ’ 


where for any fc E {1,. .., iV + 1}, (/3j) is a 2fc + 1-uple such that for all j E {1,..., 2fc + 1}, one has 

2fc+l 

0 < |A| < ••• < \fhk+i\ < |a| < N and ^ fa = a + (0,2), 

j =i 

and C(^.) is a constant. 

Assume first that |/3 2 fc+i| = N. Then for any j £ {1,..., fc}, |/L| < 2. It follows 

U>lP2fc + l/-| ^ 1 

I W 2 


1 \JM\ 


1 


^(|^CL 2 ,oo)|^ 2fc+1 C 


< 


c'dd^jici 


*B„-! - Bo 1 / 2 


Bo 1 li2 Bo 1/2 

Now, if |/? 2 fc+i| = N — 1, then either |/5 2 fc| = 3 and |/5j| = 0 for any j < 2k — 1, or \/3j\ < 2 for 
any j < 2fc. The latter case is estimated as above, while in the former case, one has 


Bo 


r (A)| < 

r k 11,2 ^ 


Bo 1 / 2 


C(\d x (\ w 2 ^\d^(\ wL Jd^d x (\ L2 < )|C|^ 


Otherwise, one has j3j\ < N — 2 for any j £ {1,... ,2k + 1}, and in that case, 


Bo ' 5 l^-Bo C(|C 


x T 


)|C 


X 1 


Altogether, this yields for N > 4 


(B.25) 


Bo 


r (a) l < 
r 5 Il2 - 


Bo 1 / 2 


C(Bo 1 ,^e 2 ,|CL. N )|C 


Finally, we notice that 


c<a) def 7 + i y 3 ga f A 

" “ Bo J x V(l + Me 2 |5xC| 2 )5/2 X V 

may be estimated, when l<|a|<_/V—1, as 

(B.26) \d x a a C {&) \ L2 <^ 2 C(Bo-\^ 2 ,\(:\ xN )|C|*m (l<\a\<N-l). 

Bo- 1 Bo" 1 

The definition of the operators a, b, c, r^ and estimate (B.8) follows from (B.lO)-(B.ll I, 


p3T2} - p3T3| |BT7t - pT8| ), p^22| l- |lL23} , ( |B^ -( |R25| as well as ( |B.26) when 1 < \a\ < N- 1. 
Estimate (B.9) is obtained identically, using in particular the trivial estimates 


|-ff( e Ci) — roc < -Hi,/i 0 e|Ci — Cs | Loa , |^( e Ci) ~ H(e ( 2 )\ L2 < H ltho e\(i - C2 


IL 2 ’ 


and 


|#«i) _ H(e( 2 )\ L1 < \H(e(i) - H(e( 2 )\ H1 < C(H 2 ^ ho , |eCi|wi.°°’ | e C21 w i,oo)e|Cl — C 2 \ h i, 
where H and H t ho are as in Lemma [B.3| This concludes the proof of Lemma B.4 □ 
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B.3 Preliminary results 


In this section, we prove that the operator a[e£, eiu] (resp. b[e£]), introduced in Lemma B.4 is 
symmetric, continuous and coercive with respect to the space Xg o _! (resp. Y p° M ), provided that 
some conditions are satisfied by (e£, eu>). These requirements can be seen as sufficient conditions 
for the hyperbolicity of the system, and permit to control the energy solutions to the quasilinear 
system for positive times (Section B.41, and eventually prove the well-posedness of our system 
(Section B.5). 


Lemma B.5. Let G Hf x Z p M be such that ef satisfies ( |B.7| ) with /i 0 > 0. Then one has 

a[eC,ew] £ £(X° o _!; (X^.J*). b[e£] € £(Y F 0 „; (Y F ° )*) and c[e£,ew] € £(Yp° ; (5^°)*). Moreover, 
there exists Kq,K\ = C(m, hff , e|C|) such that 


V/,5eX° o _ X) 
V/,se^, 
V/,set 
v/ g y f °, 


|<a[eC ,ew]f , g) {x o _ i} *| < Xi(l + T F \wf zl )|/| 


x° 


\9\xo 


!(&[<]/ , I < -K’iI/Ii^I^Iy-o, 

I<c[eC ,ew]f , 5)(yo,).| < eR\\w\ zl \f\ y0 \g\ y0 , 


1 


(&[<]/ , /)(y F 0)* > j^|/|y F 0 . 
Assume additionally that there exists kg > 0 such that 
(B.27) (7 + S) — e 2 max {(hf 3 + 7 /i( _ 3 )|u'| 2 } > k 0 > 0. 

Then there exists K,K' 0 = C^m, 1 , k^ 1 , e I Cl H3 ) such that if 


Yp tu 


< K 


-1 


V / e X Bo-i> (a[e(,ew]f , f) (x 0 j* > j^\f\ x o ■ 

B° iln Bo- 1 


(B.28) 

then 


Proof. We establish each result for /, 5 G 5 (R) so that all the terms are obviously well-defined and 
in particular the (X* — X) duality product (with X = X 3 or Yp-j,) coincides with the L 2 scalar 
product; the result for f,g G Xg o _ t or Yp° M is then obtained by density of <S(R) in X^.j and Yf ,,, 
and continuous linear extension. 

One has, after integration by parts, 

(B.29) (b[<]/, 5 ) 0 -). = (&[<]/, ff) L2 

= f W^ /ff + ^^(^F 2 {/ l 2 - 1 /})(5,F 2 {ft. 2 - 1 5 }) + dx. 

^ 1/12 o o 

It follows easily 

I (&[<]/ , g)(y° )*| < Xi|/| 0 |ff|yo ■ 

I V / I 1 p/x j pM 

We write again for the coercivity inequality, 


(b[<]f , = f ^ 

Jr /; -i 


hl + Hr l l/|2 + f /i 2i^F 2 {/ l2 - i /}i 2 + ^/ l ?ia x F 1 {/ l r l /ii 2 a*. 


It follows immediately, since e<( satisfies 

1 + 7 


Mo 


1 


m)f, f)w > r yp^zri/| L 2 + YI^ F2 ^“ i /}| L2 + ^i^F 1 {/ ir i /}| i2 - 
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Now, by Lemma [B. 2 1 one has 

\d x ^if \ L 2 = |a c F 1 {(l + eC)^r 7}| L2 < (1 + e|C| z o )\d x F f}\ 

and similarly for \d x F 2 {hf l f}^ L2 . We conclude 

\f\ 2 Y ^ < C'(m,/i- 1 ,e|C| z o M ) x (b[e(]f , f) (Y ^y- 

By similar argumentation, one easily shows that the operator c[e£, ew] is well-defined and con¬ 
tinuous from U F ° fl to (Y F(i )*, and satisfies the third estimate of the statement. 

We show now the coercivity of a[e£, ew] under additional assumption (B.271. We write 
(a[e(,ew]f , f) {x o _j* = (a[eC, ™]/. f) L * 


(7 + S) - e- 


hi + 7^2 
( h\h 2 ) 3 


H 2 I/I 


7 + S 


\d x f\ 2 


with 


Ri = 
1 


(hi(d x F!t{hT l w}) 2 - -whfd x ^{hld x Ft{h^w}}) x /,/ 

+ ^(h!d x ^{hr 2 wf},d x ^{h- 2 wf}) L2 - 2(a,F ?{(h- 2 w) x /} , (h^F^w}) x z)^. 


Bo (1 +/xe|a x C| 2 ) 3/2 


i 2 


dec + /ie 2 (i? 2 - 7 -Ri), 


Using Cauchy-Schwarz inequality and Lemmata |B.2| and |B.3| one has the following estimate 
/re 2 |i? 2 - 7-Ril < t 2 \w\zi C '( m > h o S e |C| J?3 )|/|^o 

<T f | W |^ C(m,/io 1 ,e|C| H3 )|/|io , 


where the last identity follows from Lemma [B. 1 1 
From (B.27), one has immediately 

(u[eC, ew]f , f) {x o - /ie 2 (i ?2 - 7-Ri) > min ^ h 


o, 


7 + 5 


(1 + H^Clr J 3/2 . 


x I/I 


A' 0 


The existence of Kq,K such that (B.28) implies 


|/| 2 y o < K' 0 (a[e(,ew\f , /> (A -o _ 

Bo -1 Bo 

is now straightforward. 

One shows similarly that a[e£, ew] : A'g o _! —> (A”g ,)* is well-defined and continuous, and 
satisfies the first estimate of the statement. This concludes the proof of Lemma [B. 5 [ □ 


The following Lemma is a direct consequence of Lemma |B.5| 

Lemma B.6. Let (£, w ) T £ Hf. x Z F(1 be such that satisfies (B.7). Then b[e£] : is 

a topological isomorphism with: 

V/et |(b[eC]) _1 /|ym < K 0 \f\ {Y ^) t , 

with Kq as in Lemma m 

If, additionally, (ef,ew) satisfies (B.27)-(B.28), then a[e£,eu;] : A 3 — > (A° _,)* is a topolog¬ 

ical isomorphism with: 

|(a[eC,ew])~ 1 /| A o < K 'o\f\ (x o 

Rn" 1 V Ro 1 ' 


with Kq as in Lemma B.5 
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Proof. By Lemma B.5 b[e£] : Y F j, —>• (Y F )*)* is well-defined, continuous and coercive. We then 
deduce by the operator version of Lax-Milgram theorem that b [e£] is an isomorphism from Yj? onto 
(Y FfI )*. The continuity of the inverse follows from the continuity and coercivity of b[e£]: 


|&K] f\ Y ^ < K 0 {b[eQb[eQ f , b[e(] < K 0 \f\ {Y ^ y |b[eC] f\ Y ^- 

The whole discussion is identical for a[e£, ew\, replacing Y F ° M with A'g i, and Kq with K' 0 . 
We conclude this section with the following result. 


□ 


Lemma B.7. Let (£, w;) T £ Hf. x Z Ffl be such that e( satisfies (B.TI. Then the operator a[e£, ew] : 
-+ (X° _o * is symmetric: 


"Bo” 


V/,3 G-X^bo-h (a[eC,H/ , = <£»[<> ew]g , 

The same result holds true for b[e£] and c[e£, ew\, replacing Jg o _i with Y F(J . 


Proof. The symmetry property for b[e£] is straightforwardly seen from (B.29). The other operators 
require a slight rewriting. In particular, notice 

di7 li[hi,w]» = [hi(d x Ff{h~ l w}) 2 - ^w/q“ 3 <9 x Ff {/if<9 x Ff x • 

+ (A -2 ™) x d x Ff{(h 2 d x Ff{h~ l w}) x •} - (h^F^h^w}) x d x Ff{(h~ 2 w) x •} 


- -{hfw)d x ^{h\d x Ft{{hfw).}} 


and 


(dQ^NH + d 2 n f ;[h i ,w])» = (hi r 2 <tf; x 


- x .} + (h^F^h^w}) x d x F fiK 1 X .} 

+ h- 1 x d x Ff{hld x ^{{h^w) x .}} + l -(h~ 2 w) x d x F>f{h\d x Ftihr 1 x .}} 


are obviously symmetric, since d x Ff is skew-symmetric. The result is now clear. 


□ 


B.4 A priori estimates 

We now consider the quasi-linearized system arising from Lemma |B.4| 


(B.30) 


dt( + d x u> = r 1 , 

bd t w + d x a( + d x a a C + c d x w = r 2 , 


where we denote for conciseness a = a[e£, ew\ (and similarly for a a , b, c), as defined in Lemma B.4 
and r 1 ,?" 2 are remainder terms to be precised. More accurately, we introduce a regularized version 
of (B.30). Denote J„ = (1 — vd 2 ) -1 / 2 and consider 


(B.31) 


dtC + J ld x ui = r 1 , 

bd t w + J 2 3 x aC + J 2 9 x a Q C + i v ti v d x w = r 2 . 


We obtain below a uniform a priori control of the energy of any solution, and then, in Lemma |B.10| 
a similar estimate on the difference between two solutions. 
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LemmaB.8. Lett] ^ (C, w) T , U d = (C,w) T 6 i°°([0, T}- x Y& ), C G L°°([0,T]; (X^.,) 2 ) 
and r = (r 1 ,?- 2 ) 1 ^ £ L 1 ([0,T); x (Yp^)*) satisfying (B.31) with v £ [0,1]. Assume moreover 

that U(t) satisfiesand ( |B.28[ ) with h 0 1 ,k 0 1 ,K 1 uniformly for t £ [0,T]. Then one 
has 


E°(Uy 12 < Co(S 0 (t /| t=0 ) 1/2 + (^ 2 )||C|| i 0 o ([ 0 , T ] ;(A -O o _ i) 2 ) )e At + Co f f{t')e x ^d£, 


with 


A = Co x + T F \\w\ 
,-i 7.-1 


L=°([0,T];Z2 ) 


)> /(*) 


= 7" 


A'° ,X(Y» 


+ M e c 




flTic? Cg C(m,/ 1 q ,fc 0 ,Ki IMLooqoT]-J f 4 ,xYp 4 ,,)^* 


Remark B.9. 77ie energy estimate is uniform with respect to v £ [0,1], and holds in particular 
for solutions to the non-regularized system ( B.30). 


Proof. Since U, U £ I“([0,T];lg , x Y p ,,), all the components of equation (B.31) are obviously 


r4 X Y 4 

Bo -1 X L 

well-defined in L 2 . We compute the L 2 inner product of the first equation with aC + a a C, and 
add the L 2 inner product of the second equation with w. Recalling that a, b, c are symmetric (by 
Lemma B.71, and since J„ is symmetric and d x is skew-symmetric, we obtain after straightforward 
manipulations 


( B - 32 ) ^ {\( a CX) L 2 + {Ca a c) L2 + ^(bw,w) L2 ^ = 

^([<9t,a]C,C) i2 + (C,d t (a a ()) L 2 + b]w,w) L2 + ^([d x , c] i v w, i u w) L2 

+ (r 1 , aC + o a C ) L 2 + (r 2 , w) L2 • 

We estimate below each of the components of the right-hand-side. These estimates follow from the 
product estimates of Section |B.1[ as in the proof of Lemma |B.4| For the sake of conciseness, we do 
not detail all calculations but rather provide the precise estimates for each component. 

(/) = f ([di,a]C,C) t 2 - O ne has, by definition, 


[d t , a] C = -e 2 C 9 t (G(eC)M 2 ) - ( [dt, diK^hz, w]] C - l[d u [/ii, w]] C 


1 + 5. 

"bT 




(1 + pe 2 \d x (\ 2 ) 3 / 2 


4 


where G(ef) = ^ t and 7 Z^[hi,w\ is defined in (|B. 19). 

(^1^2) - 


The first contribution is easily estimated: 

( - e 2 (d t (G(eO\w\ 2 ),() L2 < C(m,/io\ |eC| wli00 ) x 


21 |2 
e \w 




The third component is estimated after one integration by parts: 


( _ i^ ( a ‘ (u+^ci 2 r 3/2 Rc), c) 


< C(m, \C\ W2 


) x fie 2 


1 

Bo 


dxC 


2 

L 2 ’ 


Treating the last contribution is more involved, as [dt, d±7l^[hi .“'l] is the sum of many terms. 
However, all of these terms may be dealt with as in the proof of Lemma [B.4| using integration by 
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parts if necessary, one may ensure that the operator d x F( L applies only once to each ( and since 


much regularity is assumed on £ £ X'k , Lemmata 


B.2 


B.3 


and then Lemma 


B.l 


yield 


MeV l (\dt,diT^[hi,w]]CC) L2 < C 0 x eVli^JCIyo, < CoT F H|fjClxo 




with Co = C( m, h 0 , |£| 4 )• Altogether, we proved 

A Bn ~ 1 


(B.33) 


|(/)| <C(m,h^,\C\ xi )x /xe 2 + T F | W |; 2 |C|^o • 

Bo -1 ' / Bo -1 


(77) = (\d t , b]w,w) L2 . One has, by definition, 

' hi + 7/12 


[d t ,b]w = d t ( 1 Ji ^ 2 )w + ne(dQl[h 2 ](w) -'ydQ F 1 [hi\(w))dtC, 


where d Q 1 ) is defined in (B.16). The first term is estimated as 




L 2 


< C( m, h 0 , 9tC roo) x e w 


• i 2 


L 2 ’ 


For the second term we have after integration by parts and by triangular inequality 

(dQ 4 F N(w)a t C,w )^ ^ \(h?d x F(‘{hr 1 w} , d x F^{h- 2 d t (w}) L2 \ 

+ | (hUdtOd^ih^w}, d x P}{hT'w}) L2 | 

By Lemmata |B.2| and |B.3[ one immediately deduces 


(B.34) 


\(II)\<C(m,h 0 L ,\C\ X i ) x e\w\y 0 ■ 


(III) = f ([9a;, c] J „w, J vw) L2 . One may proceed similarly as above, and one obtains without any 
additional difficulty 

(B.35) |(7//)| < C(m, h^ 1 , |C| w , \w\ yi ) x e|j„ii;| 2 0 < C(m, 1 , |CL 4 ,\w\ Yi ) x e\w\ 0 . 

Bo -1 1 F^ 1 FI L Bo -1 1 F^ I FI L 


(IV) = (C, d t (a a ()) L 2 ■ After one integration by parts, one has 

(n rs o 27 + ^ (~ i a f(d x d e i()(d x ()(d x (j 

(IV) = -3pe [ dx(,d t 


Bo 


( l + ^2\ dxC \2)V2 JJ L2 ■ 


je{ 1 , 2 } v 

Recall ( = (Co>Ci) T G X^ x Xg l5 so we easily deduce by Cauchy-Schwarz inequality 
(B.36) |(717)| < C(m, |C| W 3 ,oo) x /re 2 |C| x o |CUx )2 - 

Bo " 1 V Rn 1 ' 


(V) d = (r 1 , aC + a a C) L2 + ( r 2 ,w) L2 . By Lemma 


B.5 


one obtains 


(B.37) |(R)| < C(m,/i 0 1 ,fc 0 1 ,e|C| ff3 )k 1 | :Y o |C| x o + H^C(m 1 \C\ w2t00 )\r 1 \ 0 |C|, 0 )2 

x Bo - 1 Bo - 1 Bo - 1 V Bo 1 
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Altogether, plugging (B.33),(B.34)(B.35 I,(B.361,(B.37) into (B.32) yields 
(B.38) 

^ Q(aC,C) L 2 + (C,a a C) i2 + < Co (e + T F \w\ 2 z ^) E°(U) + CoC^U) 1 ' 2 , 

with Co = C(tn, h^ 1 ,E 4 (U)), C± = |r| x0 _ i x( Y -o i )*+Me 2 |C|( X i _ i)a , and E°(U) = |C|^ 0 _ i +|“ , |yo i - 
By Lemma B.5 there exists Ko,K\ = C(m, ,kQ X ,K, E 4 (U)) such that 


(B.39) 


Let us now estimate 


JLpO/fr' 1 ' ' 1 

Kr 


E°(U) < -(aC,C ) L2 + ^{bw,w) L2 < K 1 E°{U). 


1 


(C,aaC ) L 2 < fJ-e 2 C 2 x |C|_ Y0 |C| a .o < -p,e 2 C 2 x (M + M\(\ x0 ) 


with C 2 = C(m, \d x C,\ w ^ t00 ) and arbitrary M > 0. Choosing M = jit 2 C 2 K {1 , (B.39) yields 
(B.40) 

where we denoted 


^E\U) - M < E°(U) < (Ki + -L-)E°(U) + M, 
ZKq ZKq 


M = f max ^K 0 ^e 2 C 2 ) 2 \C\ x0 ^ [>, E°(U) ^ £(oC,C)+ (C,a«C) ia + ^w,w) L ,. 


The differential inequality (B.381 may therefore be reformulated as 
d 


^(£°([7) + M) < 2K 0 C 0 (e + T f |uj| 


) (e°(U) + M) + y/2K~ 0 CoCi (E\U) 


-M 


1/2 


We deduce 


(E°(U) + M) 1/2 < (E°(U | t=0 ) + M) 1/2 e At + C 0 [ C , i(t , )e A(t_t ' ) dt', 

Jo 


where A, C 0 are as in the statement of the Lemma. Using (B.401 and augmenting C 0 if necessary, 
the energy estimate is now straightforward. □ 


Lemma B.10. Define two tu ple o f solutions to (B.30I, (Ui,U\,ri) and (U 2 ,U 2l r 2 ), satisfying the 
same properties as in Lemma B.8 (with Ci = C 2 = 0J- Then one has 

£7°(f7! - f7 2 ) 1/2 < C 0 E°(U 1 | t=0 — U 2 | t=0 ) 1/2 e A * + C 0 

Jo 


with 


A = C 0 X 


(e+T F 


\ Wl \\L°°([ 0 ,T]-,Z^) 
-1 7,-1 


). f(t) 


= ri-r 2 


XL-!><(*&) 


+e\U 2 \ 


(w£ 


\Ui-U 2 \ 


and Cq C(m,h 0 ,k 0 ,K, ||^i|| i «,([ 0)T ] ; j i: 4 _II^ 2 IIl«>([o,t];X^ _ lX v F ^))‘ 
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Proof. The difference between the two solutions satisfies the system 
( d t ( C 1 -C 2 ) + d x {w! -w 2 ) = r\-r\, 

{ bidtiwx - w 2 ) + SccditCi - C 2 ) + tid x {w\-W 2 ) = r\ - r\ + r d iff, 
where we denote cq = <x[eCi,eWi] (and similarly for b*, c*), and 


3 

rdiff d = (b 2 ^ bi)d t W 2 + (d x a 2 - c^di)^ + (c 2 - d)d x w 2 d = y~Vdiff- 

i—1 


The Lemma is a straightforward consequence of Lemma B.8 (with v = 0), once rdiff is estimated. 
We focus on the most difficult term, namely = ( d x a 2 — 9 x ai)C2- 
Let / £ Yp 0 ^. One has 


r ditfJ) L 2 = [ ((G(eC 2 )|u ; 2 | 2 - G(eCi)|^i| 2 iC2) 

J K 

- ne 2 fd x (^&i1ll[eC, 2l w 2 } - di^Ki. w i] + 7<ii'£i[ e C2, w 2 ] - ydi^ifeCi,wi])6 


I ±*fd 2 

Bo J x 


(1 + He 2 \d x ( 2 \ 2 ) 3 / 2 (l + pe 2 |^Ci| 2 ) 3 / 2 


d x C 2 


where G(e£) = 


def h\ + jhl 


(hih 2 ) 3 

Since |G(eC 2 — e Ci )| L 2 < e|C 2 — Ci | L 2 x su Pye[eC 2 ,eCi] d is straightforward that 

\dx ((G(eC 2 )k 2 | 2 - G(eCi)|wi| 2 )C 2 )| i2 < G 0 (e|Ci - ( 2(^1 + \wi - w A m ) |C 2 | w i.=o, 

with G 0 = G(m, h^ 1 , |Ci| ioo , |C 2 | i0 o, \wi\ Lao ,\w 2 \ Loo ). 

Similarly, 


Bo 


^l\^x (((1 + ^ e ~\dx(,2\ 2 ) 3,/2 — (1 + H e2 \d x (l\~) 3 ^ 2 )dxC2^j\ L 2 

< Ht 2 Co\d x C,i — d x C, 2 \ xl \d x ( 2 \ w 2,oa, 

Rn-1 x 


with G 0 = G(m, \ d xC,i\ w i,o°)- 


I w£ 


As for the last component, recall di7?. F is defined in (B.19). Proceeding as in the proof of 
Lemma IB. 41 we obtain 


7 2 Ve 2 |(di^[eC 2 ,^ 2 ]C 2 - di^, F [eCi,Wi]C 2 ,/) L2 | 

< e 2 Go x (Iwi — w 2 \ ! + e|Ci — C 2 L 1 ) IC 2 L 1 l/k-o 

1 1 p/i. 1 '^FA 4 1 1 FA 6 

with G 0 = G(m, /ig k ICiLa, IC2L3. | u ’iLi , | W2 Li )■ 

1 x ^FA 4 ^FA 4 

Altogether, we find 

I (2) I I * I I I 

rdiffl(Y^)* - eC ° I b\w *’°°\ U2 - Ul \ *| 0 _ lX y F v 

with G 0 = G(m, h^ 1 , k^ 1 , K, |tAi | 4 4 , |C/ 2 | 4 y4 )• 

Bo” 1 _ Bo -1 FAt _ 

All the other terms in rdiff are estimated in the same way, and Lemma |B.10| now directly follows 
from Lemma IB. 81 □ 
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iteration scheme. Uniqueness, and continuous dependence with respect to the initial data follow 
from Lemma TB. 101 


Lemma B.ll. Let w, (, r 1 , r 2 £ H°°([0, T] x M) be such that ( B.7 ),( B.27 1 , ( |B.28 ) hold. Then for 
any U° ((' 0 ,w 0 ) T £ iJ£°(]R.) 2 , there exists a unique solution U *= (£, ii>) T £ iL°°([0,T] x K) 2 


satisfying (B.30) and U\ t=0 = U° ■ 


Remark B.12. One could assume only continuity in time and finite (but large enough) regularity 
in space on w, C, r. 


Proof. We first consider the regularized system introduced in ( B.31| | and that we rewrite (recall 
that, by Lemma B.6 b^ 1 : (Y^,)* — > Y)), is well-defined and continuous) as 


(B.41) 


d t tu + J 2 v d x Wv = r 1 , 

dtWv + b~ l S 2 d x aC, + b _1 J 2 <9 x a a C + b -1 J„cJ v d x w v = b^ 1 r 2 . 


Since J„ '= f (1 — vd 2 ) - 1 / 2 is of order —1, (B.41 1 is a system of ordinary differential equations on 
X® _! x Yff , which is solved uniquely by Cauchy-Lipschitz theorem. More precisely, for any v > 0 
and r = (r\r 2 ) £ C°([0,r];X» . x (y°. )*), C € C°([0, T\; (Xg ,) 2 ) and U° £ X° x 7°,, 


there exists a unique U v = f (f u ,w v ) T £ C' 1 ([0, T]; Xg o _i x Yj?^), solution to (B.411 with initial data 

f7| t=0 = (C°,^°) T - 

Differentiating N times (B.311 and proceeding as in the proof of Lemma [~B.4[ one can check that 
d x U v satisfies (B.311 with obvious modifications to r 1 ,r 2 and (. Thus, by the above argument, 


d x tl v £ C' 1 ([0, T]; Xg o _, x U F ° M ), and it follows (since N may be chosen arbitrarily large) that 
U v £ C ,1 ([0,T];iL^ o (]R)). In particular, d t U v | t=0 £ Hf°. 

Applying the above argument to dtU „ after differentiating (B.311 with respect to time, one 
deduces d t U v £ C 1 ([0, T\; iJ£°(lR)), a nd b y induction U v £ H°°([0,T] x 


Applying the estimate of Lemma B.8 to d x U v with N £ N given, one has 

E°(d?U u ) < M, 


- 1 kT\K,T,E 0 (d?U 0 ),\\((,w,(,r) 
Let us now consider V,, = U,, — U„>. V,, „> satisfies 


with M = C(m, h 0 , n, 0 




), uniform with respect to v > 0. 


.31) with C = 0, V v y | t=0 = 0 and 


r v y = (J„ - j V ')d x w v 


r lv' = (4 - ji')dxaC,v' + (j^ - il,)d x a a C, + (J, y cJ„ - J^cJ v >)d x w v >. 


= 1 and J,. — J 


—> 0 (v —> i/) and thanks to the above 
one deduces that U v 


H ‘-¥H L 
0 (v —► v'). By Lemma 


B.5 


Since for any s £ K, || -U 
energy estimates, one has r„ „/ „ 0 

’ IA Bo -1 X VF 

is a Cauchy sequence of C°([0, T\; Xg o _, x Y F(1 ). Therefore there exists a limit that we denote 
U £ C°([0, T]; A’° o _ 1 x Yp^), which satisfies the non-regularized (i.e. v = 0) system, namely ( B.30| ). 

The above energy estimates on d x U v being uniform with respect to one has U £ L°°([0, T]; Hff). 
By (B.30(, we deduce d t U £ L°°([0, T]; iJ£°), and by induction U £ if°°([0,T] x R). 

Uniqueness of the solution follows when applying the energy estimate of Lemma |B.8| to the 
difference between two solutions. □ 
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We can now conclude with the proof of our main result, Theorem |5.3| 
Proof of Theorem \5.3\ Let us define Friedrichs modifiers, = 1(|£)| < x) and 

U„ U = U° = f {(C?“j 2 nC°, ^j 2 ^ 0 )}| Q |<iV- 


For each n > 1, we define, thanks to Lemma B.ll U„ {(Cn Q \ wi a ^)}| a I < tv as the unique solution 
to u « L 0 = u ° as well as ( |B.30[ ), where (using the notations and definitions of Lemma |B.4| 
a = a[e£ n _!, eto n _i] and similarly for b, c, r = r^ a ); a£ = 0 if |a| < N — 1 and af = a a QAi 

o 


otherwise. Our iteration scheme is initialized with smooth and time-constant Uo = Up. 

Lemma B.ll defines at each step U n G C([0, T n \; iJ£°), where 

T n (h' 0 , k' 0 , K\ M') d = max |t > 0, such that E N (( n , w n ) 1/2 < M'E N (C /°) 1/2 

and (£„,«;„) satisfies (B.7),(B.27),(B.28) with h' Q ,k' 0 , AT'j. 

One has T n > 0 as soon as h' 0 < ho, k' 0 < ko, K' > K and M' > 1 , by standard continuity 
arguments. Le t us prove that T n can be bounded from below, uniformly with respect to n G N. 

By Lemma B .8 we have the energy estimate for U^ *==' w^) T with any |a| < N: 

E 10 (C/to)) 1/2 < C 0 (E 0 {Ui a) \ t=0 ) 1/2 + ne 2 M'E N (U°) 1/2 )e Xt + C 0 f /(f')e A ^-*'W, 

Jo 


for any t G [0, T n -i(h' 0 , k' 0 , K', M')\ and with 
A = Cq x ( e + T f |K_ 


(e +T F ||u;„_i|| io 0 ([o iT ); Z 2 ^)), f(t) |?’ ( 




+ /ie 2 M 1 E N (U 0 ) 1 / 2 , 


and where Co = C(m, (h' 0 ) 1 , (k' 0 ) 

Notice that d a U n 7 ^ but one can check (differentiating the equations satisfied by U n ) 

that d a U n satisfies (B.30) with a remainder term [e( n , ew n , eCn-i, erOn-i] which is estimated 
identically as in Lemma B.4 This yields, for any t G [0, min{T„_i, T n }], 

E N (U n ) 1 / 2 < C 0 e xt E N (U°) 1 / 2 (l + C , 0 M'tx (e + T* /2 |+ T F |w„_ x |^)), 

with A, C 0 as above, and C' n = C(m, (Hq) -1 , M', E n (U 0 )). 

We deduce that there exists M*, ^ = C(m, h ^ 1 , k ^ 1 ,K, E N (U 0 )), independent of n, such that 


T n (h 0 /2, /co/2,2 K, M*) > T*/X', X' d = e + T/ 
and that for any t G [0, T*/A'], one has 

(B.42) E N {U n ) 1/2 < M* E n (U°) 1/2 . 

Let us now consider V n = U n — U n -\. Notice first that 

E j 04 | t=0 ) = Ei{{U n - U n _f) | t=0 ) < 2- 2 ^ n ~^E n (U°). 

One can control E°(V n ) from Lemma 


1 r Y' I 0 I Z 
1 + 7 F | IV | Z 2 , 


B.10 


using the above, the estimate on r 1 °^ — r l?-i given by 


Lemma B.4 as well as the energy estimate (B.42). Similar estimates on d a V n for 0 < |a| < 2 yield 

E 2 (V n ) l/2 < C 0 2-" (Ar - 2 ) e Vt + C 0 A' [ ^; 2 (l/ ^ _l) 1 / 2 e V(^ - t ' ) d^ , , 

Jo 
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with Co, A' as above. Therefore, restricting T$ < T* if necessary, the sequence U n = U° + J2j=i b 
converges in C°([0, T^/X']; x Y^). 


Using that U n is uniformly bounded in C°([0, T^/A']; x Ypl), the logarithmic convexity 

of Sobolev norms yields that U n converges strongly in C'°([0, T**/A']; x Yp^ -1 ). The limit 

U = lin Woo U n belongs to L oo ([0,T#/A , ];^_ 1 x Y$) n C°([0, T#/A'[; X*~\ x Y^" 1 ) and then 
by classical argument belong to (^([O, TyA']; x Y F )(). It is now straightforward to check that 

U satisfies system (B.30), and therefore (by Lemma B.4| (5.1). 

By passing to the limit the energy estimate (B.42), one deduces the energy estimate of the 
statement. The uniqueness of the solution is a consequence of Lemma |B.10| applied to the difference 
between two solutions (see also Proposition 5.8 1 . Theorem 5.3 is proved. □ 
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